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elsewhere.  .Similarly,  the  principles  involved  in  assessing  blast  effects  in  cold  regions  arc  the 
same  as  those  that  apply  to  blast  effects  generally.  The  reports  therefore  summarize  princi¬ 
ples  and  data  for  general  explosions  technology,  and  then  present  the  procedures  and  data 
for  cold  environments  within  this  framework. 

The  purpose  of  the  scries  is  to  provide  a  convenient  reference  source  for  engineers  faced 
with  pro  ilems  of  explosions  or  blasting  in  cold  regions.  Because  not  all  engineers  are  well  ac¬ 
quainted  with  explosions  technology,  relevant  physical  principles  are  explained  or  summar¬ 
ized,  but  there  is  no  attempt  to  explore  the  underlying  theory  in  depth,  nor  is  there  any  treat¬ 
ment  of  the  practical  aspects  of  explosives  use  and  .safety  practices.  These  topics  are  covered 
well  in  Army  Technical  Manuals  and  Army  iTateriel  Command  publications,  as  well  as  in 
commercial  blasters’  handbooks  and  in  textbooks. 
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Blasting  and  Blast  Effects  in  Cold  Regions 

Part  III;  Explosions  in  Ground  Materials 


MALCOLM  MELLOR 


INTRODICTION 

Considerable  attention  has  been  given  to  explo¬ 
sive  excavation  in  cold  regions,  and  to  the  effects 
produced  by  explosions  in  frozen  media  such  as 
snow,  ice  and  permafrost.  Many  experiments, 
field  tests  and  theoretical  studies  have  been  carried 
out,  but  the  results  that  have  been  produced  are 
insufficient  to  give  a  complete  and  coherent  pic¬ 
ture  if  frozen  materials  arc  treated  in  isolation.  To 
gain  a  full  understanding  of  the  relevant  behavior, 
the  limited  data  for  frozen  materials  have  to  be  in¬ 
terpreted  in  the  context  of  general  explosion  tech¬ 
nology,  as  it  applies  to  unfrozen  rocks  and  soils. 
In  this  report,  explosion  effects  in  solid  ground 
materials  are  summarized  in  general  form,  with 
frozen  materials  treated  as  special  types  of  soils 
and  rocks. 

Ice  is  a  characteristic  constituent  of  frozen 
ground  materials,  and  it  has  physical  properties 
that  differ  considerably  from  those  of  typical 
rock-forming  minerals  and  soil  grains.  For  exam¬ 
ple,  ice  density  is  about  one-third  that  of  silicates, 
and  the  compressibility  of  ice  under  extreme  pres¬ 
sure  i.i  much  higher  than  the  compressibility  of 
typical  rocks.  However,  in  spite  of  the  unusual 
properties  of  ice,  frozen  materials  such  as  snow, 
ice  and  frozen  soil  behave  much  like  other  soils 
and  rocks  when  they  are  subjected  to  blasting  or  to 
blast  effects. 

For  most  practical  purposes,  ice-bonded  soils 
can  be  treated  like  rock  or  concrete  when  design¬ 
ing  explosive  excavation  techniques  or  when  as¬ 
sessing  blast  effects.  For  example,  water-saturated 
frozen  sand  is  like  sandstone,  frozen  silt  is  like  silt- 
stone,  and  water-saturated  frozen  gravel  is  like 
concrete.  Massive  ice  also  behaves  like  rock  under 
explosive  attack.  It  attenuates  a  stress  wave  more 
strongly  than  do  typical  hard  rocks,  but  its  crater¬ 
ing  characteristics  are  not  much  different  from 
those  of  common  sedimentary  rocks  like  lime¬ 
stone:  the  idea  that  ice  is  a  fragile  material  that 
shatters  easily  is  quite  unrealistic.  Snow  deposits 


have  low  bulk  density  and  high  porosity  (com¬ 
monly  40%  to  80%),  and  it  is  no  surprise  to  find 
that  snow  is  highly  effective  in  attenuating  stress 
waves.  However,  the  cratering  characteristics  of 
dense  snow'  are  not  much  different  from  those  of 
common  rocks  and  soils,  since  low  strength  and 
low  density  are  to  some  extent  offset  by  strong  at¬ 
tenuation. 

When  the  pore  water  freezes  in  common  rocks 
and  minerals,  there  is  a  significant  increase  in 
strength.  This  can  affect  the  economics  of  a  large- 
scale  mining  operation  that  is  efficient  and  well- 
optimized,  but  it  is  not  of  major  significance  in  or¬ 
dinary  blasting  operations  for  excavation  and  con¬ 
struction. 

This  report  is  concerned  mainly  with  explosions 
inside  the  ground  material.  The  effects  of  explo¬ 
sions  at  or  above  the  ground  surface  are  consid¬ 
ered  in  more  detail  in  Part  I. 


DEEP  UNDERGROUND  EXPLOSIONS 

When  a  concentrated  explosive  charge  is  set 
deep  underground  with  the  borehole  stemmed,  the 
initial  effects  of  the  explosion  are  spherically  sym¬ 
metrical  as  long  as  the  surrounding  material  is  not 
strongly  anisotropic  (Fig.  1).  The  ground  is  stressed 
by  an  initial  shock  wave  and  also  by  expanding 
gas.  There  is  local  heating  by  adiabatic  compres¬ 
sion  of  the  ground  material,  and  the  surface  of  the 
blast  cavity  is  exposed  to  the  high  temperature  of 
the  explosion.  The  effects  on  the  surroundings 
vary  to  some  extent  with  the  properties  of  the 
ground  material.  For  example,  highly  porous  ma¬ 
terial  with  air-filled  pores  can  compact  locally  to  a 
high  degree,  whereas  very  dense  material  can  suf¬ 
fer  gross  displacement  only  by  distributing  strains 
and  displacements  over  a  large  surrounding  vol¬ 
ume. 

When  the  shothole  is  completely  filled  by  the 
charge,  the  ground  material  that  is  in  direct  con¬ 
tact  with  the  charge  experiences  very  high  pres- 
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Figure  1.  Effects  of  u  fully  contained  explosion  on  the  surrounding 
solid  material.  Spherical  symmetry  for  a  spherical  charge:  cylindri¬ 
cal  symmetry  for  a  long  cylindrical  charge. 


siires  (approaching  detonation  pressure).  Since  the 
initial  amplitude  of  the  shock  wave  from  a  typical 
high  explosive  far  exceeds  the  deviatoric  yield 
strength  ol  any  solid  material,  the  ground  material 
close  to  the  charge  undergoes  intense  compression 
that  is  essentially  hydrodynamic  and  adiabatic. 
Br  ;tle  Triatfinl  such  as  rock  is  completely  pulver¬ 
ized  in  this  zone  (Fig.  1).  As  distance  from  the 
charge  increases,  the  shock  wave  amplitude  de¬ 
creases  due  to  geometrical  spreading  of  the  wave 
front  and  also  due  to  dissipation  in  the  ground 
material.  Inelastic  (plastic)  compression  of  the 
material  becomes  progres.sively  less  severe,  and  the 
shear  resistance  becomes  increasingly  important. 
At  sufficient  radius  from  the  charge,  the  behavior 
of  the  material  is  that  of  an  elastic  solid  rather 
than  the  hydrodynamic  behavior  of  pulverized 
material.  In  the  transition  zone,  the  rock  exhibits 
an  elastic-brittle  response,  and  radial  cracks  are 
formed  by  tensile  hoop  stresses,  either  from  the 
radially  propagating  shock  or  from  the  pressure 
pulse  produced  by  gas  pressure  in  the  explosion 
cavity.  In  principle,  cracks  can  be  pressurized  and 
propagated  Py  the  expanding  gases,  provided  that 
the  gases  are  not  “sealed  in”  by  intense  pulveriza¬ 
tion  (and  possibly  fusion)  of  the  cavity  wall.  How¬ 
ever.  crack  propagation  by  gas  intrusion  probably 


does  not  occur  in  a  fully  contained  explosion.  Be¬ 
yond  the  zone  of  cracking  the  stress  wave  propa¬ 
gates  elastically  at  the  acoustic  velocity  for  the  ma¬ 
terial,  producing  dilatation  and  shear. 


SIZE  OF  THE  CHAMBER  AND 
THE  CRUSHED  ZONE 
(Spring  Hole  Blasting) 

The  cavity  produced  by  a  deeply  buried  explo¬ 
sion  is  termed  variously  a  chamber,  a  camouflet,  a 
spring  hole,  or  simply  a  cavity.  “Chamber”  is  a 
geneiai  term,  “caM«.  /‘let”  is  an  old  military  term 
for  an  explosion  smothered  by  burial,  and  “spring 
hole”  usually  means  a  cavity  blasted  at  the  base  of 
a  shothole  to  provide  space  for  more  explosive. 
The  limit  of  the  cavity  is  not  always  easy  to  define, 
since  fragmented  material  can  collapse  back.  Nor 
is  the  full  extent  of  the  crushed  zone  easy  to  meas¬ 
ure,  since  it  is  extremely  expensive  to  cut  cross  sec¬ 
tions  through  hard  rock  at  great  depth. 

Test  data  for  rocks  and  soils 

In  various  test  programs  carried  out  by  the  U.S. 
Bureau  of  Mines,  spring  holes  were  blasted  by  con¬ 
centrated  charges  and  loose  fragments  were  then 
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Figure  2.  Sununary  of  spring  hole  data  from  USBM  tests  in  granite. 
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Figure  3.  Results  of  USBM  spring  hole  tests  in  salt.  (After 
NichoUs  and  Hooker  1962a.) 


blown  out  through  the  shothole  by  compressed  air. 
These  studies  produced  some  of  the  best  North 
American  data  for  small  explosions,  but  it  has  to 
be  recognized  that  the  volume  of  material  excavat¬ 
ed  is  somewhat  arbitrary.  The  volume  measured  is 
certainly  more  than  the  volume  of  open  cavity 
when  no  fragments  are  removed,  but  it  is  likely  to 
be  significantly  less  than  the  total  volume  of  rock 
that  is  subject  to  severe  cracking. 

Figure  2  summarizes  the  results  of  four  series  of 
tests  made  in  granite  by  the  U.S.  Bureau  of  Mines, 
using  a  variety  of  explosives.  This  plot  of  volume 
against  charge  weight  displays  considerable  scat¬ 
ter,  but  there  is  a  systematic  linear  trend  over  al¬ 
most  5  orders  of  magnitude.  The  indication  is 
that,  for  typical  explosives  well-coupled  to  gran¬ 
ite,  the  volume  of  thoroughly  crushed  rock  is 
about  0.14  ftVlb,  or  8.7  x  10  ’  mVkg.  There  is  un¬ 
certainty  by  a  factor  of  2  on  either  side  of  this 
value  (the  volume  could  be  twice  as  much,  or  half 
as  much).  The  volumes  plotted  represent  the  total 
volume  of  the  excavated  cavity,  with  no  deduction 
for  the  original  volume  of  the  charge. 

Figure  3  gives  the  results  of  another  series  of 
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tests  made  in  salt.  Here  the  range  of  the  data  is 
smaller,  hut  the  data  points  lie  very  close  to  a 
straight  line,  with  little  scatter.  The  yield  in  salt 
was  only  half  that  in  granite— about  0.064  ft'/lb, 
or  4  '  10  '  m '/kg. 

One  other  study  in  rock  produced  a  good  linear 
correlation  between  chamber  volume  and  charge 
size.  In  Project  Buchanan  (Gillespie  1972),  spring 
charges  were  fired  in  a  weak,  weathered  granite, 
giving  the  results  shown  in  Figure  4.  The  data  that 
have  been  plotted  here  are  for  spring  charges  in  a 
fresh  hole;  during  the  project,  second  and  third 
charges  were  loaded  successively  into  spring  holes 
produced  by  previous  detonations.  The  data  in 
F  igure  4  give  specific  cavity  volume  as  approxi¬ 
mately  0.06  ft’/lb.  or  3.75  v  l()  '  m'/kg. 

Other  test  programs  in  tonalite,  weathered  clay 
shale,  sandy  overburden,  and  clay  have  produced 
data  that  do  not  fit  a  linear  relation  between  vol¬ 
ume  and  charge  weight.  In  general,  these  tests 
cover  relatively  small  ranges  of  charge  size,  and 
they  cannot  be  taken  as  strong  evidence  against 
the  linear  relation,  especially  as  the  log-log  plots 
suggest  power  relations  with  exponents  both 
above  and  below  unity.  Using  mid-range  values  to 
establish  representative  yields  for  these  tests,  spe¬ 
cific  cavity  volumes  are  about  0.06  ft'/lb  (j.Vx 
10  ’  m'  kg)  for  tonalite,  2.2  ft'/lb  (0.14  m'/kg) 
for  weathered  clay  siiale,  I  .k  ft'Tb  (0. 1  m'/kg)  for 
clay  ,  and  1.2  It'  lb  (0.07  m'  kg)  for  sandy  over¬ 
burden. 

Hcnrych  (19^4)  accepts  without  question  that 


chamber  volume  should  be  proportional  to  charge 
weight  for  concentrated  explosions  inside  rocks 
and  soils.  For  fully  contained  nuclear  explosions 
in  rock,  Glasstone  and  Dolan  (1977)  also  accept 
approximate  proportionality  between  chamber 
volume  and  the  energy  yield  of  the  explosion. 
From  the  evidence  summarized  above  there  seems 
to  be  no  good  reason  to  doubt  this  proportional¬ 
ity,  so  that  chamber  size  can  be  expressed  as  a  spe¬ 
cific  volume  (Table  la,  lb). 

If  it  is  assumed  that  the  cavity  produced  by  a 
compact  charge  is  approximately  spherical,  the 
volume  measurements  give  an  effective  radius  for 
the  cavity.  The  specific  volume  data  for  granite  in 
Figure  2  imply  a  mean  scaled  radius  of  0.32 
ft/lb  ',  or  0.13  m/kg  '.  This  is  just  over  two 
charge  radii — 2.1,  2.2  and  2.3  charge  radii  for 
charge  specific  gravities  of  1.0,  1.2  and  1.4,  re¬ 
spectively.  The  results  for  salt  imply  a  scaled  radi¬ 
us  of  0.25  ft/lb  '  or  0. 10  m/kg  '.  For  charge  speci¬ 
fic  gravity  of  I.O,  1.2  and  1.4,  the  cavity  radius  is 
1.6,  1.7  and  1 .8  charge  radii  respectively.  The  data 
for  weak,  weathered  granite  suggest  a  cavity  radi¬ 
us  of  0.24  ft/lb  '  or  0.10  m/kg  '.  This  is  equiva¬ 
lent  to  1.6,  1.65  and  1.7  charge  radii  for  charge 
specific  gravities  of  1.0,  1.2  and  1.4  respectively. 
In  tonalite  (Fig.  5),  cavity  radius  is  about  0.24 
ft/lb  '  or  0.1  m/kg  ',  i.e.  the  same  as  weathered 
granite.  These  results  are  summarized  in  Table  la. 

Henrych  (1979)  gives  cavity  radii  for  hard  rocks 
in  the  range  1 . 1  to  2.5  charge  radii,  or  0.06  to  0.13 
m/kg  (Table  Ic).  For  weak  and  weathered  rocks. 
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valuer  arc  in  the  range  1.7  to  4.0  charge  radii,  or 
0.04  to  0.21  m  '-.g  '. 

In  soils  and  very  weak  roeks,  the  cavities  are 
much  bigger.  In  weathered  clay  shale  the  radius  is 
about  0.81  I't'lb  or  about  0.32  m/kg  ’.  For  the 
slurry  explosive  that  was  used  (DBA  22M),  this  is 
ec|uivalent  to  6  charge  radii.  In  saturated  clay,  cav¬ 
ity  radius  is  roughly  the  same,  and  in  sandy  over¬ 
burden  a  little  less  (Table  la).  Henrych  (1979) 
gives  data  for  a  variety  of  soils.  The  overall  range 
of  cavity  radius  is  3.8  to  13.1  charge  radii,  or  0.2 
to  0.7  m  kg  ’  (Table  Ic). 

Glasstone  and  Dolan  (1977)  give  representative 
values  of  C7'  ity  radius  for  fully  contained  nuclear 
explosions  deep  inside  rock  (>2000  ft,  or  >600 
m).  For  dense  silicate  rocks  such  as  granite,  the 
cavity  radius  is  35  ft  kt  (11  m/kt  ')•  For  dense 
carbonate  rocks  (limestone,  dolomite)  the  cavity 
radius  is  25  ft/  kt  '  (7.6  m/kt  ')■  If  these  values  are 
converted  directly  into  conventional  HE  units,  as¬ 
suming  I  kt  fully  equivalent  to  1000  tons  of  TNT, 
the  radii  are  approximately  0.27  ft/lb' '  (0.11 
m/kg  ')  and  0.2  ft/lb  '  (0.08  m/kg  ')  for  silicate 
and  carbonate  rocks  respectively  (Table  lb).  If  a 
cratering  equivalency  factor  is  used  in  the  conver¬ 
sion,  say  1  kt  =  lO”  lb  TNT,  the  cavity  radii  be¬ 
come  0.34  ft/lb  '  (0.14  m/kg  ')  and  0.25  ft/lb  ' 
(0.098  m/kg  ')  for  silicate  and  carbonate  rocks  re¬ 
spectively.  These  values  are  in  good  agreement 
with  the  HE  values  in  Table  la. 

In  theoretical  studies  of  rock  blasting,  Kutter 
and  Fairhurst  (1971)  postulated  an  “equivalent 
cavity"  of  fractured  material  formed  by  the  shock 
wave,  and  pressurized  by  the  expanding  gas.  They 
deduced  that  this  equivalent  cavity  should  have  a 
radius  of  6  charge  hole  radii,  or  6  charge  radii,  for 


spherical  symmetry.  l  or  cylindrical  charges,  the 
corresponding  figure  was  9  charge  hole  radii,  or  9 
charge  radii.  The  postulated  equivalent  cavity  is 
therefore  bigger  than  actual  cavities  I'ormed  in 
hard  rock. 

Test  data  for  cylindrical  cavities  blasted  by  lin¬ 
ear,  or  columnar,  charges  are  not  easy  to  come  by. 
Table  Id  gives  a  few  reported  values  for  cylindri¬ 
cal  cavities  in  soils.  There  are  also  some  reported 
values  of  about  3  to  4  charge  radii  for  the  “crushed 
zone"  when  cylindrical  charges  are  fired  in  hard 
rock.  These  can  probably  be  interpreted  as  cavity 
radii. 

Experimental  values  for  specific  chamber  vol¬ 
ume  can  be  interpreted  as  powder  factors  (i.e.  re¬ 
ciprocal  of  specific  hole  volume,  which  is  equal  to 
specific  charge).  The  values  for  chambering  in 
rocks  are  more  than  two  orders  of  magnitude  high¬ 
er  than  the  powder  factors  that  apply  to  typical 
blasting  operations,  but  the  values  for  chambering 
in  soils  and  very  weak  rocks  are  quite  low  .  For 
rocks,  the  effective  powder  factor  is  about  200- 
450  Ib/yd'  (115-270  kg/m’)  and  for  very  weak 
rock,  soils,  and  frozen  soils  the  range  is  about 
12-50  Ib/yd’  (7-30  kg/m’). 

Spring  holes  in  frozen  materials 

Spring  hole  measurements  have  been  made  in 
frozen  soils,  usually  without  any  attempt  at  exca¬ 
vation  of  loose  debris.  Benert  (1961)  fired  60% 
gelatin  dynamite  in  a  frozen  till  described  as 
“gravelly  silt.’’  The  results  (Fig.  6)  give  a  general 
idea  of  spring  hole  volume  in  frozen  till,  but  they 
cannot  be  used  to  define  a  relation  between  vol¬ 
ume  and  charge  weight,  since  the  length/diameter 
ratio  of  the  charges  .ncreased  with  charge  weight, 
reaching  a  value  of  12  for  the  biggest  charges. 

Other  tests  were  made  in  frozen  silt,  using  sensi¬ 
tized  nitromethane  and  composition  C-4  (Mellor 
1972,  Sellmann,  unpublished).  These  were  not  sys¬ 
tematic  studies,  and  the  results  (Fig.  6)  are  not 
suitable  for  defining  general  relations  between 
chamber  volume  and  charge  weight.  Bauer  (1975) 
fired  compact  charges  of  ANFO  (S.G.  0.85)  in 
frozen  clay.  The  results  follow  a  linear  trend,  al¬ 
though  Bauer  himself  did  not  accept  a  linear  rela¬ 
tion  between  cavity  volume  and  charge  weight. 

Looking  at  the  overall  compilation  of  data  in 
Figure  6,  it  appears  that  spring  hole  volumes  for 
compact  charges  in  frozen  silt  and  frozen  clay  may 
be  quite  similar,  with  specific  volume  about  0.53 
ftVlb  (3.3  X  10  ■'  m'/kg).  The  equivalent  radius  for 
a  spherical  chamber  is  about  0.5  ft/lb  ‘  (0.2 
m/kg  ').  This  is  3.2  charge  radii  for  a  charge  speci- 
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Figure  6.  Summary  of  spring  hole  data  for  frozen  Figure  7.  Summary  of  stress  attenuation  data  for  spherical 
soils.  propagation  in  granite. 


fic  gravity  of  1 .0,  or  3.6  charge  radii  for  a  specific 
gravity  of  1 .4, 

The  data  for  frozen  till  are  probably  representa¬ 
tive  of  the  performance  of  a  compact  spring 
charge  at  the  lower  end  of  the  data  range,  but  at 
the  upper  end  of  the  range  they  represent  cylindri¬ 
cal  rather  than  spherical  geometry.  In  this  materi¬ 
al,  specific  volume  is  perhaps  about  three  times 
smaller  than  in  frozen  silt  and  frozen  clay. 


SIZE  OF  THE  CRACKED  ZONE 

It  is  not  easy  to  determine  the  extent  of  radial 
cracking  around  a  deeply  buried  explosion.  To  get 
a  feel  for  the  size  of  the  cracked  zone,  results  of 
the  few  direct  observations  have  to  be  combined 
with  deductions  from  indirect  evidence. 


Research  studies  (Olson  et  al.  1973)  suggest  that 
radial  cracks  from  a  concentrated  charge  in  gran¬ 
ite  may  extend  about  1 8  to  20  charge  radii  (2.5-2.75 
ft/lb'^’,  or  1. 0-1.1  m/kg'''’).  Another  estimate  can 
be  made  by  considering  the  empirical  stress  atten¬ 
uation  curve  for  granite  in  Figure  7.  If  the  com¬ 
pressive  radial  stress*6T  the  graph  is  assumed  equal 
to  the  tensile  circumferential  stress,  and  if  the  ten¬ 
sile  strength  of  the  granite  is  taken  as  2000  Ibf/in.' 
(13.8  MPa),  then  radial  cracks  could  perhaps  ex¬ 
tend  out  to  a  radius  of  5.5  ft/lb (2.2  m/kg'  ’),  or 
about  39  charge  radii.  If  the  same  type  of  argu¬ 
ment  is  applied  to  the  attenuation  of  peak  radial 
strain  (Fig.  12),  and  appropriate  values  are  taken 
for  the  tensile  failure  strain  of  the  rock,  another 
estimate  is  obtained.  Taking  the  tensile  failure 
strains  for  granite  and  limestone  as  4.5  x  10  "  and 
3  X 10  "  respectively,  the  indicated  values  of  maxi- 
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mum  crack  radius  are  about  29  charge  radii  for 
granite  (4  ft/lb' or  1.6  m/kg'"’)  and  about  40 
charge  radii  for  limestone  (5.5  ft/lb'"'*,  or  2.2 
m/kg' ').  A  very  crude  upper  limit  estimate  can  be 
made  by  drawing  on  attenuation  curves  for  peak 
particle  velocity,  together  with  empirical  damage 
criteria  for  near-surface  waves.  A  safe  limit  of 
particle  velocity  for  internal  cracking  of  sound 
rock  is  about  25  in./s  (635  mm/s);  for  a  spherical 
wave  in  granite,  this  level  is  reached  at  a  radius  of 
about  13  ft/lb  '  (5.2  m/kg  '^'),  or  about  94  charge 
radii. 

These  rough  estimates  are  summarized  in  Table 
2a.  It  seems  fairly  clear  that  maximum  cracking 
radius  is  in  the  range  10  to  100  charge  radii.  It  is 
not  unlikely  that  the  actual  range  is  about  20  to  40 
charge  radii  for  spherical  symmetry. 

For  /cing  cylindrical  charges,  there  are  more  di¬ 
rect  observations  on  the  size  of  the  cracking  zone 
(Cattermole  and  Hansen  1962,  Siskind  et  al,  1973, 
Siskind  and  Fumanti  1974).  Maximum  radius  of 
the  cracked  zone  seems  to  be  in  the  range  15-55 
charge  radii.  This  fits  quite  well  with  empirically 
established  values  of  hole  spacing  in  practical 
bench  blasting.  A  typical  spacing/diameter  ratio 
of  45  for  bench  blasting  implies  a  crack  extent  of 


about  45  charge  radii.  The  extreme  limits  in  bench 
blasting  are  20  to  80.  Looking  at  established  prac¬ 
tice  for  smooth  blasting  and  presplit  blasting,  and 
assuming  a  decoupling  ratio  of  3.5,  the  ranges  of 
hole  spacing  in  charge  diameters  are  45  to  56  for 
smooth  blasting  and  28  to  42  for  presplitting.  This 
can  perhaps  be  interpreted  as  an  expectation  that 
crack  extension  from  a  single  charge  will  be  in  the 
range  of  45-56  and  28-42  for  smooth  blasting  and 
presplitting  respectively. 

Another  estimate  can  be  based  on  the  static 
plane  strain  elastie  analysis  for  a  hole  in  a  thick 
plate.  Consider  a  hole  of  radius  t?  in  a  plate  that 
has  in-pJanc  principal  stresses  a,  =  o,  =  P.  If  the 
hole  has  internal  pressure  p,  the  radial  and  circum¬ 
ferential  stresses  at  radius  r  from  the  center  of  the 
hole  are,  respectively; 


Of  =  P-i-  {p-  P){a/rY 

(1) 

o„  =  P  -  [p-  P){a/ry. 

(2) 

Up  »  P, 

o„  =  -  Or  =  -p(a/ry. 

(3) 

8 


When  a  borehole  is  pressurized  abruptly,  an 
elastic  stress  field  can  establish  itself  in  the  sur¬ 
rounding  medium  at  the  acoustic  velocity  for  the 
medium.  For  that  part  of  the  material  that  does 
not  fracture  during  passage  of  the  shock,  the  limit¬ 
ing  radius  for  crack  formation  is  the  value  of  r 
where  equals  the  tensile  strength  of  the  rock  Oj. 
This  is  assuming  that  the  borehole  pressure  can  be 
sustained  for  longer  than  it  takes  to  establish  the 
elastic  stress  field. 

A  first  inclination  is  to  identify  hole  radius  a 
with  the  drill  hole  radius  d/2,  d  being  the  drill  hole 
(and  charge)  diameter.  Corespondingly,  the  hole 
pressure  p  might  be  identified  with  the  detonation 
pressure  of  the  explosive  or  with  some  fraction 
of  Pp.  However,  it  may  be  more  reasonable  to  take 
a  as  half  the  diameter  of  the  blast  cavity,  or  half 
the  diameter  of  the  crushing  zone,  either  of  which 
can  be  denoted  Kd,  where  A’  >  1 .  The  correspond¬ 
ing  hole  pressure  p  would  then  be  the  pressure 
which  is  reached  after  adiabatic  expansion  of  the 
explosion  products  from  the  detonation  pressure 
Pq  at  diameter  d  to  p^^  at  diameter  Kd. 

For  the  first  set  of  assumptions,  the  potential 
limit  for  crack  extent  is  given  by 

r  =  aip^/a^f'  =  (d/2)(p^/aj)''\  (4) 

For  the  second  set  of  assumptions,  the  limiting  ra¬ 
dius  for  crack  extent  is  given  by 

r  =  {Kd/2)(p^^/aj)''^  (5) 

where,  for  a  long  cylindrical  charge,  adiabatic  ex¬ 
pansion  gives 

Ph=Pd^~^'' 

in  which  7  =  1.2  for  expanding  explosion  pro¬ 
ducts.  Thus 

r  =  (Kd/2)(p^/K^^ 

=  K'-'(d/2)(Po/ajy^\  (7) 

With  A  =  2  to  3,  A‘"'  is  in  the  range  0.8  to  0.87, 
so  the  choice  of  assumptions  does  not  have  a 
strong  effect  on  the  final  result. 

For  exploratory  estimates,  take  the  tensile 
strength  of  fairly  strong  rocks  as  1000  to  2000 
lbf/in.‘  (6.9  to  14  MPa).  Take  A  =  2  to  3,  and  p„ 
in  the  range  0.7  x  lO*  to  2.8  x  lO'  Ibf/in.’  (5  to  19 
GPa).  The  factor  (Pi/a^)  ^  controls  the  result; 
pairing  high  and  low  limits  for  Pp  and  Oy  gives  a 


range  from  19  to  53,  while  the  pairing  of  low  and 
high  limits  for  both  reduces  the  range  to  26  to  37. 
These  are  maximum  crack  lengths  in  charge  radii 
for  the  first  and  simplest  set  of  assumptions.  For 
the  second  set  of  assumptions,  maximum  crack 
lengths  are  80%  to  87%  of  these  values. 

Table  2b  summarizes  the  deduced  limits  for  max¬ 
imum  crack  length  around  a  long  cylindrical 
charge  in  rock.  Again  it  seems  safe  to  conclude 
that  maximum  crack  length  will  be  in  the  range  10- 
100  charge  radii.  The  actual  range  is  probably  nar¬ 
rower,  say  20  to  55  charge  radii. 


ATTENUATION  RELATIONS  FOR 
SHOCK  PROPAGATION 

Spherical  symmetry 

In  reviewing  the  available  data  for  spherical 
propagation  of  blast  effects  in  ground  materials, 
the  standard  forms  of  presentation  used  in  the  lit¬ 
erature  will  be  followed  (next  section).  In  the  stan¬ 
dard  presentations,  data  are  plotted  on  logarith¬ 
mic  scales  for  good  reason,  since  the  values  can 
span  many  decades.  There  is  also  a  common  prac¬ 
tice  of  representing  the  data  by  straight  lines  on 
log-log  plots,  and  hence  by  simple  inverse  power 
relations,  i.e.  I  =  h  r'",  where  /  is  the  magnitude 
of  the  effect,  r  is  the  radius  from  the  source,  and  /» 
and  n  are  empirical  constants.  This  is  convenient, 
but  it  is  not  necessarily  good  science;  n  has  to  vary 
for  different  ranges  of  r,  /  must  tend  to  infinity  as 
r  tends  to  zero,  and  there  is  no  separation  of  geo¬ 
metrical  effects  and  dissipation  effects. 

As  a  stress  wave  propagates  soherically,  it  atten¬ 
uates  because  the  wave  spreads  over  an  ever-in¬ 
creasing  surface  area,  and  also  because  energy  is 
dissipated  in  displacing  and  straining  the  medium 
through  which  the  wave  passes.  The  effect  of  geo¬ 
metrical  spreading  is  the  same  for  any  material 
and  at  any  distance  from  the  source;  the  area  of  a 
spherical  wave  front  is  proportional  to  By  con¬ 
trast,  the  effects  of  dissipation  tend  to  decrease 
with  distance;  close  to  the  source,  much  energy  is 
expended  in  pulverizing  and  fracturing  rock,  but 
at  distant  range  the  wave  is  elastic  and  not  subject 
to  much  loss. 

For  a  simple  wave,  the  energy  is  proportional  to 
the  square  of  the  amplitude.  The  area  traversed  by 
the  wave  front  is  a  spherical  surface,  and  thus  the 
energy  is  spread  over  an  area  that  increases  with 
r\  The  energy  per  unit  area  is  then  inversely  pro¬ 
portional  to  r\  The  amplitude  is  thus  inversely 
proportional  to  r.  If  the  wave  propagation  were 


perfectly  lossless,  with  no  dissipation  of  energy  in 
the  medium  and  no  change  of  wavelength,  the  am¬ 
plitude  /  would  attenuate  as 

I  =  I.C  r,)  '  (8) 

where  /.  and  r.  are  reference  values  of  amplitude 
and  radius  respectively,  and  a  finite  lower  limit  is 
set  for  r.  In  air,  pressure  wave  attenuation  far 
from  the  source  is  not  much  more  rapid  than  \/r, 
and  in  water  the  attenuation  follows  an  approxi¬ 
mate  1/r  trend  even  at  fairly  close  range. 

In  general,  it  is  to  be  expected  that  energy  will 
be  absorbed,  or  dissipated,  as  a  wave  travels 
through  any  real  material.  Whatever  the  term  used 
for  the  dissipation  (internal  friction,  absorption, 
etc.),  the  expectation  is  that  the  amplitude  of  a 
plane  wave  (no  spreading)  will  attenuate  exponen¬ 
tially,  i.e. 

/  =  /o  e  (9) 

where  /« is  the  amplitude  at  r  =  0  and  cr  is  an  atten¬ 
uation  constant  with  the  dimension  of  (length)"'. 

Combining  the  effects  of  geometrical  spreading 
and  internal  dissipation,  a  rational  form  for  the 
attenuation  relation  of  a  non-dispersive  wave  in 
isotropic  material  is 

/  = /„(/•//•„)-' e""'"  (10) 

where  ro  is  a  lower  limit  of  r,  approximately  equal 
to  the  charge  radius,  and  /o  is  the  value  of  /  at  r  = 
ro. 

Equations  similar  to  eq  10  are  used  in  geophys¬ 
ics  to  express  the  attenuation  of  seismic  waves, 
and  applicability  to  explosive  stress  waves  in  rock 
has  been  demonstrated  by  Duvall  and  Petkof 
(1959)  and  by  .Atchison  and  Roth  (1961).  Equation 
10  can  be  written  in  logarithmic  form  as 


and  therefore  a  semi-log  plot  of  the  product  /r  (log 
scale)  against  r  (linear  .scale)  should  be  a  straight 
line,  with  the  slope  giving  «.  Duvall  and  Petkof 
(1959)  plotted  data  for  the  peak  amplitude  of 
strain  waves,  and  obtained  linear  trends  for  many 
different  explosives,  four  different  rock  types,  and 
various  charge  sizes.  Fogelson  et  al.  (1959)  did 
similar  work  with  six  different  explosives  in  gran¬ 
ite.  Atchison  and  Roth  (1961)  obtained  similar  re¬ 
sults  for  strain  waves  in  marble.  Values  of  the  at¬ 


tenuation  constant  a  were  found  from  the  slope  of 
the  regression  line.  All  data  were  for  the  elastic  re¬ 
gion,  starting  just  outside  the  limit  of  major  crack¬ 
ing.  Scaled  values  of  «  by  Duvall  and  Petkof  var¬ 
ied  slightly  with  the  explosive  type,  but  mean  val¬ 
ues  were: 

Granite  0.034  (ft/lb  ')  ' 

Sandstone  0.048  (ft/lb 

Chalk  0.027  (ft/lb^''’)  ' 

Shale  0.026  (ft/lb '  ')-' 

A  value  of  0.085  (ft/lb  "*)"'  was  accepted  for  marl- 
stone,  but  the  data  show  a  distinct  break  in  slope 
on  the  semi-log  plot.  Fogelson  et  al.  (1959)  found 
a  mean  value  of  0.031  (ft/lb '  ')  '  for  granite.  The 
scaled  values  listed  above  are  applicable  when  r  is 
used  in  corresponding  scaled  form.  If  r  is  ex¬ 
pressed  in  feet,  then  the  values  ot  a  given  above 
are  numerically  equal  to  values  of  a  in  ff.  Atchi¬ 
son  and  Roth  (1961)  give  a  mean  value  for  marble 
of  a  =  0.011  ft  '.  Taking  a  with  dimensions 
(length)"',  the  SI  equivalents  of  these  values  are; 

Granite  0.110  m"'  (0.103  m"') 

Sandstone  0.158  m"' 

Chalk  0.088  m"' 

Shale  0.086  m"' 

Marble  0.037  m"' 

While  eq  10  has  some  attractive  features,  its  ap¬ 
plicability  is  likely  to  be  restricted  to  a  certain 
range  of  radius,  in  the  same  way  as  the  power  rela¬ 
tion  is  restricted,  since  the  dissipation  mechanisms 
are  very  different  at  different  ranges  (e.g.  crush¬ 
ing,  fracture,  anelastic  vibration).  Attenuation  by 
internal  dissipation  is  also  dependent  on  fre¬ 
quency;  in  general,  attenuation  rate  increases  with 
increase  of  frequency.  Fligh  frequency  compo¬ 
nents  of  a  complex  wave  tend  to  be  filtered  out 
first,  and  the  frequency  characteristics  are  thus 
changed  at  long  range.  This  effect  may  be  fairly 
insignificant  in  some  materials,  but  not  in  others. 
Finally,  the  geometric  attenuation  is  only  a  1/r 
function  for  a  non-dispersive  wave,  or  a  pulse  with 
constant  phase  duration. 

Taking  a  purely  empirical  approach,  measured 
relations  between  /  and  r  can  be  described  very 
well  by  a  polynomial,  e.g. 

I  =  aj  '  +  a-r  ■  +  u.r  '  (12) 

in  which  the  coefficients  Un  are  curve-fitting  pa¬ 
rameters  and  the  radius  r  can  be  expressed  in  any 
convenient  .scaled  form.  However,  there  is  little 
point  in  fitting  a  polynomial  to  badly  scattered 


10 


data  over  a  small  range  of  r  while  ignoring  boun¬ 
dary  conditions.  Some  empirical  attenuations  ob¬ 
tained  in  this  way  give  ridiculous  values  of  /  when 
the  charge  radius  is  substituted  for  r  (e.g.  Henrych 
1979). 

W  hatever  the  relative  merits  of  different  attenu¬ 
ation  relations  for  wave  amplitude,  the  simple  in¬ 
verse  power  relation  is  still  the  most  favored  form 
in  e.xplosion  technology.  It  embraces  all  types  of 
attenuation  arising  from  geometry,  absorption, 
dispersion  and  frequency  filtering.  It  is  simple, 
convenient  and,  when  appropriately  restricted, 
quite  acceptable.  To  avoid  problems  with  dimen¬ 
sions  and  boundary  conditions,  the  power  relation 
can  be  written  as 

/  =  l,(r/r„Y''  (13) 

where  /„  is  the  amplitude  at  the  lower  limit  of  ap¬ 
plicability  for  distance,  ro.  For  lossless  propaga¬ 
tion,  such  as  might  occur  with  an  acoustic  wave  in 
an  ideal  elastic  medium,  n  =  \.  Where  there  is  dis¬ 
sipation,  rt  >  1;  in  rocks,  «  =  2  is  common  and  in 
soils  n  can  be  as  much  as  3  to  4  at  close  range. 

Cylindrical  symmetry 

With  ideal  cylindrical  propagation,  the  energy 
per  unit  area  of  wave  from  is  proportional  to  the 
radius  r  and  therefore  the  wave  amplitude  is  in¬ 
versely  proportional  to  r'^’.  This  is  the  same  as  the 
stress/distance  relation  for  a  static  elastic  stress 
field  (eq  3).  If  attenuation  is  expressed  in  the  form 
of  eq  13,  perfectly  lossless  propagation  is  repre¬ 
sented  by  n  -  Vi.  When  there  is  dissipation,  n  > 
'A  but  it  should  be  a  bit  less  than  the  correspond¬ 
ing  value  for  spherical  propagation.  Henrych 
(1979)  gives  a  value  of  «  =  1.44  for  sandy  loam  at 
fairly  long  range  (r/ro  >  30). 


ATTENUATION  OF  STRESS,  STRAIN 
AND  GROUND  MOTION 

Just  as  in  water  and  air,  the  amplitude  of  the 
spherically  propagating  shock  front  decreases  with 
increasing  radius  from  the  charge.  However,  the 
stress  field  created  in  a  solid  by  a  shock  wave  or  an 
elastic  wave  is  not  isotropic,  as  it  is  in  water  and 
air;  the  radial  and  circumferential  stress  compo¬ 
nents  are  not  necessarily  equal  outside  the  close 
range  “hydrodynamic”  zone.  In  a  fluid  there  is 
only  a  compression  wave,  whereas  in  a  solid  there 
is,  in  the  elastic  response  range,  a  compression 
wave  (P-wave)  oscillating  in  the  radial  direction. 


plus  a  shear  wave  (S-wave)  oscillating  transversely 
to  the  radial  direction.  The  front  of  a  shock  wave 
can  be  characterized  by  the  peak  radial  stress,  and 
this  can  be  plotted  against  scaled  radius  to  give  at¬ 
tenuation  curves  for  various  types  of  ground  ma¬ 
terials.  Measurements  of  shock  pressure  in  solid 
materials  are  more  difficult  and  uncertain  than 
corresponding  measurements  in  water  and  air.  In 
addition  to  the  usual  problems  of  frequency  re¬ 
sponse  and  calibration  for  the  gauges,  there  are 
problems  of  coupling,  impedance  matching,  and 
directional  resolution.  Furthermore,  measure¬ 
ments  may  be  affected  by  proximity  to  the  sur¬ 
face,  and  also  by  anisotropy  and  inhomogeneity  in 
the  ground  material. 

Attenuation  data  for  peak  radial  stress 

Figure  7  gives  a  compilation  of  results  for  the 
attenuation  of  peak  radial  stress  in  granite  when 
the  test  conditions  are  such  as  to  approximate 
spherical  symmetry.  Data  for  extremely  high  stress 
levels  refer  to  nuclear  explosions,  which  can  in¬ 
duce  close-range  pressures  greater  than  the  deton¬ 
ation  pressure  of  high  explosives.  The  nuclear  data 
are  from  tests  designated  Hardhat  (5  kilotons,  290 
m  depth).  Shoal  (12.5  kilotons,  366  m  depth),  and 
Piledriver  (50  kilotons,  460  m  depth).  Different 
authors  interpret  the  same  measurements  in  slight¬ 
ly  different  ways.  One  of  the  low  stress  curves 
seems  erroneous,  perhaps  because  of  incorrect 
conversion  of  units.  The  high  explosive  data  by 
Noren  (DuPont)  are  for  water  pressure  in  bore¬ 
holes,  which  is  not  necessarily  the  same  thing  as 
rock  stress.  However,  if  the  analysis  of  Khristof- 
orov  and  Romashov  (1967)  is  applied  to  Noren’s 
data,  the  agreement  with  the  other  results  is  excel¬ 
lent.  Taking  a  representative  line  through  all  of 
the  credible  data,  the  decay  of  stress  with  distance 
is  very  close  to  an  inverse  square  relation  (actual 
exponent  -1.92)  for  distances  out  to  about  25  ft/ 
lb'  '  (10  m/kg' '). 

Figure  8  gives  some  attenuation  data  for  spheri¬ 
cal  propagation  in  volcanic  tuff.  In  this  material, 
stress  levels  are  somewhat  lower  than  they  would 
be  in  granite  at  the  same  range.  Overall,  the  decay 
of  peak  stress  with  distance  is  very  close  to  an  in¬ 
verse  square  relation  for  distances  out  to  about  25 
ft/lb (10  m/kg'  ').  Comparable  data  for  alluvium 
are  shown  in  Figure  9. 

Smith  (1966)  gave  decay  exponents  for  a  variety 
of  materials.  Not  all  of  the  values  given  are  accept¬ 
ed  by  this  writer,  but  results  for  dry  sand,  clay,  silt 
and  loess  indicate  close-range  attenuation  with  an 
exponent  of -4  (Fig.  11). 
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Fi.i’iire  8.  Stress  attenuation  data  for  spherical  propa¬ 
gation  in  volcanic  tuff. 

Henrych  (1979)  shows  how  the  stress  attenua¬ 
tion  in  sand  varies  with  the  water  content  and  the 
air  content  (or  dry  bulk  density).  For  saturated 
sand,  the  attenuation  exponent  is  -1.05  with  zero 
air  content,  -1.5  with  5x10“  air  content,  -2.0 
with  10  ’  air  content,  and  -2.5  with  4x10'^  air 
content.  In  other  words,  saturated  sand  with  no 
air  gives  attenuation  like  water,  while  saturated 
sand  with  plenty  of  air  bubbles  gives  attenuation 
like  rock.  For  unsaturated  sand  of  dry  bulk  densi¬ 
ty  Q,  the  attenuation  exponent  was  -2.8  for  q  = 
1.6- 1.7  Mg/m’,  -3.0  for  q  =  1.52-1.6  Mg/m’, 
and  -3.5  for  q  =  1.45-1.5  Mg/m’. 

The  attenuation  of  stress  in  glacier  ice  and  in  de¬ 
posited  snow  is  illustrated  in  Figure  10.  In  ice, 
stress  appears  to  decrease  with  distance  raised  to 
the  power  -2.3,  a  decay  rate  that  is  a  bit  more 
rapid  than  rates  found  in  typical  rocks,  and  much 
faster  than  the  decay  rate  for  liquid  water  (see  Part 
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Figure  9.  Stress  attenuation  data  for  spherical  propa¬ 
gation  in  alluvium. 

II,  Fig.  1).  Within  the  range  of  available  data, 
overpressures  in  ice  are  well  below  the  values  for 
corresponding  distances  in  water.  By  comparison 
with  granite  (Fig.  11),  overpressures  in  ice  are 
lower  than  the  values  in  granite  at  close  range  ( <  1 
m/kg'''*),  but  comparable  at  longer  distances  (e.g. 

10  m/kg'"’^).  Attenuation  in  dense  snow  appears  to 
be  strong  in  the  near  field;  the  decay  exponent  is 
-3.75,  or  close  to  the  approximate  value  of  -4  that 
was  given  above  for  near-range  attenuation  in  soils. 
Within  the  range  of  existing  data,  the  peak  stresses 
in  snow  are  lower  than  corresponding  values  in 
air,  water,  ice,  rocks  and  typical  soils  (Fig.  11). 

Attenuation  data  for  peak  strain 

In  some  attenuation  studies  the  emphasis  has 
been  on  strain  rather  than  stress.  Figure  12  repre¬ 
sents  a  compilation  of  data  for  the  attenuation  of 
peak  radial  strain  with  distance  for  various  types 
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Figure  10.  Stress  attenuation  for  spherical  propagation  in  gla¬ 
cier  ice  and  deposited  snow.  (Ice  data  from  Ingram  and  Hal- 
per  I960,  Livingston  I960,  Joachim  1967.  Snow  data  from  In¬ 
gram  a"d  Strange  1958,  Ingram  and  Halper  I960,  Joachim 
1967,  Livingston  1968.) 


Figure  12.  Attenuation  of  peak  radial  strain  for 
spherical  propagation  in  granite.  (Data  from  Du¬ 
vall  and  Atchison  1957,  Atchison  and  Tournay 
1959,  Duvall  and  Petkof  1959,  Atchison  and  Pug- 
liese  1964,  NichoUs  and  Hooker  1965.) 
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Figure  / /.  Comparison  of  stress  attenuation  for  spherical 
propagation  in  water,  air,  rocks,  soils,  snow  and  ice. 
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/■/i;  Aitenuation  of  peak  radial  strain  for 

spherical  propagation  in  salt.  (Data  from  Nicholls 
1962,  Nicholls  and  Hooker  1962a.) 


of  high  explosive  in  granite.*  Beyond  the  range 
where  the  peak  compressive  strain  is  comparable 
to  the  failure  strain  of  the  rock,  the  response 
should  be  approximately  elastic,  and  related  to  the 
peak  stress  components  through  the  elastic  mod¬ 
uli.  In  this  elastic  range,  peak  strain  appears  to  be 
inversely  proportional  to  the  square  of  distance 
from  the  charge,  which  is  in  agreement  with  the 
stress  attenuation.  The  test  results  may  be  influ¬ 
enced  by  proximity  of  the  wave  to  the  surface. 

Figure  13  gives  an  impression  of  the  attenuation 
of  peak  radial  strain  in  salt.  The  original  reports 
have  the  detailed  data  points  fitted  by  a  line  that 
changes  slope  on  the  log-log  plot,  with  slightly 
steeper  slope  at  the  shortest  distances.  Overall,  the 
combined  results  follow  a  trend  in  which  peak 
strain  is  approximately  proportional  to  radius 
raised  to  the  power  -Vj. 

Figure  14  gives  strain  attenuation  data  fora  var¬ 
iety  of  explosives  detonated  in  a  limestone  (density 
2.36  Mg/m’).  The  original  authors  (Atchison  and 


•  Data  from  Nicholls  and  Hooker  (1962b)  have  been  omitted. 
They  do  not  fit  the  general  trend,  and  the  same  test  results  seem 
to  have  been  used  to  give  different  values  of  strain  in  a  later  re¬ 
port  (Nicholls  and  Hooker  196',' 


Q.  10 


Limestone 

(density  2.36  Mg/m^) 


a 

1 


lii.i 

100 


1000  (ft/lb 


10 

Scaled  Radius 


ICO  (m/kg 


Figure  ]4.  Attenuation  of  peak  radial  strain  for 
spherical  propagation  in  limestone.  (Data  from 
Atchison  and  Pugliese  1964b.) 
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Figure  15.  .Attenuation  of  peak  radial  strain  for 
spherical  propagation  in  marble.  (Data  from 
Atchison  and  Roth  1961.) 
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Figure  16.  Comparison  of  attenuation  curves 
for  peak  radial  strain  in  four  different  rocks. 
(Data  from  Duvall  and  Atchison  1957.) 


attenuation  rate  near  the  source,  seen  also  in  Fig¬ 
ure  12,  could  be  interpreted  as  an  indication  of  a 
limiting  compressive  strain  at  close  range. 

Attenuation  of  displacement, 
velocity  and  acceleration 

For  some  purposes,  the  internal  stresses  and 
strains  of  the  ground  material  are  of  less  interest 
than  the  displacement  and  its  time  derivatives,  i.e. 
velocity  and  acceleration.  This  is  particularly  true 
in  the  outer  zone  of  elastic  disturbance,  where  the 
effects  are  of  a  seismic  nature. 

If  a  wave  travels  through  a  medium,  passage  of 
the  disturbance  causes  each  particle  of  the  medium 
to  be  displaced  a  distance  .v  in  time  r,  where  t  is  the 
time  for  first  arrival  of  the  wave  at  that  particle. 
The  instantaneous  particle  velocity  is  A.x/At,  or 
d.x/dt;  the  instantaneous  acceleration  is  d'.x/dt'. 
The  particle  actually  oscillates  about  its  original 
position,  and  if  the  disturbance  is  idealized  as  a 
simple  sine  wave,  the  displacement  .v  can  be  ex¬ 
pressed  as  a  function  of  time: 

.V  =  A  sina't  =  A  sin(27r//)  (14) 


Pugliese  1964b)  fitted  power  relations  to  their  data 
so  as  to  give  an  exponent  of  -2.64  for  all  explosive 
types,  but  the  inverse  cube  trend  shown  in  Figure 
14  is  an  adequate  approximation  for  the  overall 
data  band. 

Figure  15  indicates  the  strain  levels  recorded  in 
marble  for  a  variety  of  explosives.  The  data  band 
has  been  drawn  here  to  represent  inverse  square 
decay. 

In  Figures  12-15,  the  range  of  variation  for  the 
strain  at  a  given  scaled  distance  probably  repre¬ 
sents  variation  of  explosive  type.  The  differences 
in  detonation  velocity  (or  pressure)  produce  dif¬ 
ferences  in  impedance  matching  between  different 
explosives  and  the  rock.  In  the  original  studies  by 
the  U.S.  Bureau  of  Mines,  the  charges  filled  the 
diameter  of  the  shothole  at  the  required  depth,  but 
the  explosives  were  chosen  so  as  to  cover  a  wide 
range  of  detonation  velocity,  and  hence  a  wide 
range  of  characteristic  impedance  (product  of  den¬ 
sity  and  velocity). 

Figure  16  compares  average  attenuation  curves 
for  four  different  rock  types.  Each  curve  repre¬ 
sents  an  average  for  several  types  of  explosives, 
with  detonation  velocities  ranging  from  6500  to 
21,000  a/s  (2000  to  6400  m/s).  The  left  side  of  the 
plot  represents  the  response  of  the  rock  where 
strains  are  greater  than,  or  comparable  to,  the  fail¬ 
ure  strains  for  uniaxial  stress.  The  trend  to  a  low 


where  A  is  the  amplitude  of  the  sinusoidal  dis¬ 
placement,  a'  the  angular  frequency  and  / the  fre¬ 
quency  in  cycles  per  unit  time.  The  particle  veloc¬ 
ity  V  is  thus 


V 


-j-  =  /I  a- cos  u,'/ 
dt 

=  A  Itrf  cos{2-n ft) 


and  the  acceleration  a  is 


(15) 
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sinu't 


=  -zl(27r/)’ sin(27r/f).  (16) 

The  peak  values  of  displacement,  velocity  and  ac- 


celeration  are 

■^max  ~ 

(17) 

^  ma.x  ~  ^  Try/I 

(18) 

^max  ~  ^ 

(19) 

and  the  maximum  value  of  velocity  is  ;r/2  out  of 
phase  with  the  maximum  values  of  displacement 
and  acceleration. 

Recalling  the  shock  equations,  outlined  previ- 
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ously  in  connection  with  air  blast  (Part  I),  it  can  be  ir. 
seen  that  there  is  a  relation  between  particle  veloc¬ 
ity  and  overpressure.  Thus,  in  principle,  measure¬ 
ments  of  peak  velocity  can  give  corresponding 
values  of  peak  pressure.  However,  to  make  this 
conversion,  the  relation  between  pressure  and  den¬ 
sity  for  the  ground  material  must  be  known.  This 
relation  is  the  Rankine-Hugoniot  characteristic. 

C  lose  to  the  source  of  an  explosion,  the  radiat¬ 
ing  stress  wave  is  certainly  not  sinusoidal — for  the 
positive  phase,  the  rise  of  pressure  with  time  is  ^ 
very  abrupt,  while  the  decay  with  time  is  quite  ^ 

gradual  and  approximately  exponential.  As  the  ® 

wave  propagates  outward,  the  front  becomes  less 
steep,  the  positive  phase  duration  may  increase, 
and  a  plot  of  displacement  against  time  for  the 
positive  phase  may  take  on  the  approximate  shape 
of  half  a  sine  wave.  Thus  the  sinusoidal  approxi¬ 
mation  is  not  too  bad  for  the  elastic  range  of  be¬ 
havior  beyond  the  hydrodynamic  zone. 

Ground  motion  from  deep  explosions  is  scaled 
on  the  assumption  of  spherical  symmetry,  which  is 
fully  justified  for  uniform  isotropic  material.  As 
usual,  radius  r  is  scaled  with  respect  to  the  charge 
radius  R,.,  or  to  the  cube  root  of  charge  weight 
H  '.  Particle  velocity  v  is  not  scaled,  just  as  peak 
overpressure  is  not  scaled  in  attenuation  relations. 
Because  time  intervals  for  explosion  effects  scale 
in  proportion  to  charge  radius  or  the  cube  root  of 
charge  weight  (see  Part  I),  the  accelerations  and 
displacements,  being  time  derivatives  and  integrals 
of  velocity  respectively,  have  to  be  scaled  with  re¬ 
spect  to  /?(.  or  W  Anticipating  a  linear  trend  of 
data  on  log-log  plots,  the  general  form  of  the  at¬ 
tenuation  relations  for  spherically  symmetric 
ground  motion  would  be 


where  K  is  a  coefficient  and  n  is  an  exponent  in 
each  relation. 


Vi  '/-.i 

-kq  ' )  (q  -lb 


Scoled  Rodius 

Figure  17.  Atlenualion  of  particle  acceleration  for. spheri¬ 
cal  propagation  in  granite.  iData  from  Nicholls  and 
Hooker  1965.) 

Data  for  approximate  spherical  propagation 
have  been  published  by  Nicholls  and  Hooker 
(1965),  who  fired  high  explosives  in  granite.  The 
original  data  have  been  re-scaled  with  respect  to 
charge  weight  and  plotted  in  Figures  17-19.  These 
results  refer  to  a  range  of  scaled  distance  that  is 
further  from  the  source  than  the  range  for  the 
stress  and  strain  plots.  Plots  by  the  original  auth¬ 
ors  omitted  some  data  points  and  fitted  trend  lines 
of  shallower  slope.  The  acceleration  data  (Fig.  17) 
show  a  trend  that  is  close  to  inverse  rube  decay, 
i.e.  n,  in  eq  22  is  approximately  3.  The  data  for 
velocity  and  displacement  (Fig.  18  and  19)  appear 
to  attenuate  less  rapidly,  with  n,  and  rtj  in  eq  20 
and  21  approximately  2.5. 

Henrych  (1979)  considers  that  internal  dissipa¬ 
tion  in  soils  is  negligible,  so  that  velocity  attenu- 
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Figure  IS.  .Anenuation  of  panicle  velocity  for  spherical 
propa.eaiion  in  granite.  (Data  from  .\icholls  and  Hooker 
1965.) 


ates  with  n,  =  2.  He  illustrates  test  data  that  show 
n.  =  2  for  both  soils  and  hard  rocks,  which  is 
rather  surprising  (range  was  1.5  to  150  m/kg  '). 

Attenuation  of  ground  motion  is  dealt  with 
again  in  a  later  section,  where  the  concern  is  the 
more  common  one  of  ground  motion  near  the  sur¬ 
face. 


STEMMING  OF  SHOTHOLES 

Traditional  stemming  materials  include  sand, 
fine  gravel,  crushed  stone,  or  mixtures  of  fine¬ 
grained  soil  (silt,  clay)  with  coarse  granular  mater¬ 
ial  (sand,  gravel).  The  particle  size  should  not  ex¬ 
ceed  10®^’o  of  the  hole  diameter.  Layers  of  different 
material  are  sometimes  placed  in  the  column,  e.g. 
clay  above  the  charge,  followed  by  sandy  gravel, 
followed  in  turn  by  coarse  gravel.  Special  stem¬ 
ming  material  packaged  in  bags  or  cartridges  can 
be  purchased,  and  special  bags  are  available  for 
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Figure  19.  Attenuation  of  particle  displacement  for 
spherical  propagation  in  granite.  (Data  from  S'ichoUs 
and  Hooker  1965. J 


packaging  standard  material  on  site.  Drill  cuttings 
are  also  used  to  backfill  shotholes,  especially  when 
the  holes  are  of  relatively  large  diameter  and  the 
cuttings  are  coarse.  In  underground  mining,  car¬ 
tridges  of  non-combustible  foams,  pastes  or  gels 
are  sometimes  used.  In  recent  years,  it  has  been  ac¬ 
cepted  that  water  alone  is  often  as  good  as,  or  bet¬ 
ter  than,  traditional  solid  stemming  in  both  verti¬ 
cal  and  horizontal  shotholes  (Anon.  1970,  Gilles¬ 
pie  1972,  Knudson  et  al.  1972,  Rooke  et  al.  1974). 
Water  has  appreciable  practical  advantages  over 
solid  materials,  which  have  to  be  procured,  deliv¬ 
ered  to  the  job  site,  and  loaded  into  the  hole  rela¬ 
tively  slowly.  Provided  the  hole  is  reasonably 
watertight  and  has  a  downward  inclination,  then 
direct  filling  with  water  is  cheap  and  easy.  Apart 
from  the  fact  that  some  holes  are  not  watertight, 
the  only  major  disadvantage  with  this  procedure  is 
incompatibility  with  explosives  that  have  poor 
water  resistance  (e.g.  unpackaged  ANFO).  These 
objections  to  water  stemming  can  be  overcome  by 
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packaeine  ihc  vsaU’i  iii  ^p^.^■ial  'Olt’--.ealin!i  plastic 
lubes,  ampoules  or  bays,  which  can  also  be  placed 
in  hori/onial  lu'les.  W  ater  siemmine  has  become 
important,  or  dominatit,  (or  shot-firing  in  West 
I  uiopean  emd  mines  in  recent  decades  (.-Xnon. 
Id'i)).  Spcci.il  vsaicr aillcd  ampoules  are  packed 
tightb  in  the  shothole  unh  positire  restraint;  the 
adsantages  include  suppression  of  dust,  fumes,  ig¬ 
nition  and  elected  tragments.  In  frozen  ground  or 
ttiassoe  ice  where  (he  temperature  is  well  below 
0  (  .  a  solid  plug  ot  stemming  can  be  formed  b\ 
tree/ing  water,  water  ice  shirrs,  or  wet  soils. 

1  he  length  r't  stemming  needed  to  contain  a 
cratcritig  esplosimi  drtes  ttoi  scent  to  base  been  in- 
scsiigated  s\ stematicidK  Practictil  guidelines  for 
crater  blastitig  ma\  recommetid  backfilling  to  one- 
ihird  or  to  otie  htil!  td  the  hole  depth,  which  is  not 
completeK  rattotial.  A  more  reasonable  e\pect;t- 
tioti  lor  this  situ.iiion  ts  tha>  the  required  length 
would  be  a  timction  ol  the  hole  diameter,  the 
properties  cd  the  stemming  tnttierial,  and  the 
charge  depth,  l  or  well-placed  granitlar  material 
that  tetids  to  lock  utider  pressure,  an  absolute 
krwer  limit  fi'r  the  required  length  of  stemming 
tnight  he  about  10  hole  diameters,  but  it  is  proba- 
bK  prudent  to  fill  to  at  least  20  hole  diameters.  In 
bench  blasting,  where  stemming  extends  from  the 
top  (d  the  charge  to  the  surface,  the  required 
length  of  stemming  ("colhir  distatice”)  is  i>pic;illy 
taken  as  proportional  to  the  burden,  which  itself  is 
it  multiple  td  the  hole  diameter.  These  recommen¬ 
dations  translate  to  14  tt'  28  limes  the  hole  diam¬ 
eter,  (Tnirney  et  al.  (1985)  ttnind  that  a  ratio  (d  at 
least  26  was  needed  for  complete  "bridging"  and 
settling,  \khen  using  water  for  stemming,  the  best 
practical  guidance  is  to  fill  the  hole  compleiely. 

Small  "blocking  charges"  are  sometimes  placed 
between  sections  of  inert  stemming.  1  hese  fire  at 
the  same  lime  as  the  main  charge,  pressurizing  the 
cohimn  of  stemming  so  as  to  resist  blowout. 

ZP  C'  C' 


SHAI.I.OW  SI  BSl  RFACK  EXPLOSIONS 

The  foregoing  is  concerned  mainly  with  deep  ex¬ 
plosions,  w  here  the  effects  are  fully  contained  and 
the  disturbance  spreads  in  a  symmetrical  pattern 
that  is  either  spherical  or  cylindrical,  depending  on 
the  charge  geometry,  A  different  situation  prevails 
when  the  charge  is  at  shallow  depth;  the  ground 
surface  then  breaks  out  and  releases  detonation 
products,  forming  a  crater  (Fig.  20),  a  trench,  or  a 
planar  surface,  depending  on  the  geometry  and 
layout  of  the  charge  or  charges.  Even  when  the 
charge  is  below  the  critical  depth,  w  here  breakout 
ceases,  the  ground  disturbance  (stress,  strain,  dis¬ 
placement,  acceleration)  is  affected  by  proximity 
to  the  surface  and  by  the  lack  of  three-dimensional 
symmetry  in  the  field  of  propagation. 

C  RATER  FORMATION 

For  purposes  of  discussion  and  analysis,  a 
cratering  charge  is  assumed  to  be  a  concentrated 
charge  of  high  explosive,  ideally  spherical  but  usu¬ 
ally  cylindrical.  A  cylindrical  charge  with  length 
diameter  ratio  less  than  four  produces  results  that 
are  indistinguishable  from  those  of  a  spherical 
charge,  and  length  diameter  ratios  up  to  six  are 
usually  considered  to  give  a  "concentrated" 
charge.  When  cratering  data  are  given,  it  is  also 
assumed  that  the  charge  lies  at  the  base  of  a  nar¬ 
row  drill  hole,  with  ICKl'^o  geometric  coupling,  and 
the  hole  is  assumed  to  be  adequately  stemmed. 
The  ground  surface  is  assumed  to  be  horizontal, 
so  that  gravity  body  forces  in  the  ground  and  grav¬ 
ity  forces  on  displaced  fragments  act  in  a  radially 
symmetric  pattern. 

Crater  geometry  varies  considerably,  depending 
mainly  on  the  scaled  charge  depth  and  the  proper¬ 
ties  of  the  ground  material.  However.  Figure  20 
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'Cr\t's  to  illustrate  the  getieral  features  of  a  typical 
crater.  Charge  depth  is  nieasi'red  from  the  original 
surface  to  the  center  of  gravity  of  the  charge. 
.■\fter  detonation,  some  material,  mainly  in  a  coni¬ 
cal  /one  above  the  charge,  is  completely  fractured 
and  subjected  to  large  displacements.  Fragments 
are  throun  upward  and  outward;  in  general,  some 
of  this  material  is  ejected  and  some  falls  back  into 
the  crater.  The  groutid  slightly  more  distant  from 
the  charge  is  fractured  or  otherwise  strongly  dis¬ 
turbed,  but  it  is  not  detached  from  the  parent 
mass,  ,-\t  still  greater  distance  the  e.xplosion  may 
produce  permanent  strain  in  the  ground  material, 
but  there  is  no  obvious  damage.  Beyond  this  zone, 
the  ground  experiences  transient  stress,  strain  and 
di>placement,  but  there  are  no  permanent  changes. 

■\  typical  crater  in  rock  (Fig.  20)  is  seen  as  a 
bowl-shaped  depression  with  a  rim  of  debris.  This 
open  hole  is  called  the  apparent  crater,  and  its 
sides  and  base  are  formed  mainly  by  shattered  de¬ 
bris.  If  all  the  loose  debris  is  dug  out  and  cleared 
away,  the  true  crater  appears.  It  is  bounded  by 
material  that  has  been  damaged,  but  not  grossly 
displaced.  The  true  crater  has  a  slightly  raised  rim 
where  the  original  ground  surface  has  been  heaved 
or  peeled  back.  Further  e.xcavation  or  core  drilling 
reveals  a  transition  zone  in  which  the  damage  or 
permanent  strain  becomes  progre.ssively  less  severe 
with  increasing  distance  from  the  surface  of  the 
true  crater.  Some  military  texts  use  very  elaborate 
schemes  for  defining  zones  of  disturbance  beyond 
the  boundary  of  the  true  crater,  and  they  also  em¬ 
ploy  a  highly  detailed  set  of  characteri.stic  crater 
dimensions.  For  most  practical  purposes  other 
than  research,  these  intricacies  can  be  ignored; 
there  are  no  sharply  defined  boundaries  between 
zones  that  experience  different  degrees  of  strain. 
Figure  20  gives  a  representation  that  is  adequate 
for  most  purposes. 

The  process  of  crater  formation  starts  with 
transmission  of  a  shock  wave  from  the  charge  to 
the  surrounding  material,  as  in  a  deep  under¬ 
ground  explosion.  The  efficiency  of  shock  trans¬ 
mission  from  the  charge  to  the  ground  depends  on 
the  geometric  coupling  and  the  impedance  match, 
as  described  earlier  for  deep  explosions.  With  effi¬ 
cient  shock  transmission,  the  ground  material  im¬ 
mediately  adjacent  to  the  charge  undergoes  in¬ 
tense  compression  and  it  is  completely  pulverized. 
This  material  is  also  displaced  radially,  forming  a 
cavity  that  is  considerably  bigger  than  the  original 
charge.  This  cavity  is  pressurized  by  expanding  gas 
from  the  explosion,  .Around  the  crushed  zone, 
cracks  radiate  outward,  some  of  them  reaching  the 


ground  surface,  or  coming  close  to  it.  As  the 
shock  radiates  outward  it  reaches  the  surface  and 
reflects. 

The  shock  is  the  first  disturbance  to  reach  and 
displace  the  ground  surface.  The  retlection  of  the 
shock  as  a  tensile  pulse  causes  small  particles  and 
larger  spall  fragments  to  break  loose  and  jump  in¬ 
to  the  air.  The  vertical  ejection  velocity  decreases 
with  the  horizontal  radius  from  surface  zero,  as 
described  previously  for  the  surface  effects  of 
underwater  explosions  (Part  II).  If  the  explosion  is 
fairly  deep  (by  cratering  standards),  this  surface 
spalling  only  affects  a  relatively  thin  surface  layer. 

Radial  cracks  probably  reach  the  vicinitv  of  the 
ground  surface  later  than  the  shock;  if  they  form 
by  extension  from  the  pressurized  blast  cavity, 
their  propagation  speed  is  much  lower  than  the 
shock  velocity  (about  20%  of  the  acoustic  P-wave 
velocity).  These  cracks,  and  their  radial  limits  rela¬ 
tive  to  the  ground  surface,  help  to  define  the  limits 
for  the  final  conical  breakout  of  fragmented  ma¬ 
terial. 

While  the  shock  produces  only  transient  stress 
waves  of  limited  thickness,  the  pressurized  blast 
cavity  establishes  behind  the  intial  shock  an  over¬ 
all  stress  field  that  is  more  akin  to  a  static  elastic- 
stress  field.  This  stress  field  is  influenced  by  prox¬ 
imity  to  the  surface;  there  is  yielding  towards  the 
surface,  and  the  blast  cavity  expands  preferential¬ 
ly  in  the  direction  of  the  surface. 

Growth  of  the  gas  bubble  displaces  and  sepa¬ 
rates  the  cracked  material;  the  surface  rises  in  a 
mound  and  accelerates,  with  velocities  greatest 
near  the  center,  thus  producing  progressive  steep¬ 
ening  of  the  mound.  As  the  ground  domes  up¬ 
ward,  the  large  tensile  strains  in  the  convex  area 
tend  to  form  two  orthogonal  crack  systems,  one  in 
the  radial  (meridional)  direction  and  the  othe-  in 
the  circumferential  (latitudinal)  direction.  These 
cracks  increase  in  width  progressively  as  the  sur¬ 
face  layer  pushes  up  and  stretches.  The  gas  then 
vents  through  the  fissures  of  the  rising  mound, 
typically  forming  a  fountain  of  faster  ejection 
plumes  that  shoot  out  of  the  top  of  the  dome. 
Small  fragments  carried  in  these  plumes  may  be 
fired  to  great  height  or  over  long  horizontal  dis¬ 
tance.  The  bulk  of  the  material  usually  drops  back 
into  the  crater  and  onto  the  area  immediately  sur¬ 
rounding  the  crater. 

■A  charge  lying  on  the  ground  surface  does  not 
make  a  big  crater,  but  it  can  fire  debris  a  long  way  . 
Burial  of  the  charge  at  shallow  depth  produces  a 
large  increase  in  crater  size.  By  increasing  the 
charge  depth,  the  surface  radius  of  the  crater  e\cn- 
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iuall\  lL■achc^  a  inaxiimini.  1  his  is  a  Ivpool  siak-i 
likcK  lo  be  ai.l\ aiuaeaiuis  in  military  eimiiK’criii”. 
\i  sumevsliai  erealei  charge  tioplli,  the  volume  o| 
the  Hue  crtitei  fetiches  ti  tnavimum,  without  much 
rciluctioii  ol  the  suilace  rtidius.  I  liis  gives  maxi¬ 
mum  elTicietiey  wlieit  hrokeii  material  is  exeavaieil 
bv  a  secoudary  mech.atiical  process.  Still  luriher 
incretise  ot  charee  depth  subdues  the  eruption  to 
the  point  vvhere  fragmented  tnaierial  heaves  up  in¬ 
to  the  tiir  atul  theti  Hops  back  into  the  crater,  leav¬ 
ing  a  moutid  ol  bulked  material  on  top.  I  bis  type 
of  crater  is  likely  to  be  advantageous  in  civil  cott- 
siruciioii,  since  it  is  deep  and  steep,  and  there  is 
very  little  flyrock.  At  slightly  greater  chtirge 
depth,  crater  formation  ceases.  This  is  termed  crit¬ 
ical  chariic  dep/h.  The  charge  forms  a  subsurface 
cavity,  as  described  earlier,  and  in  competent  ma- 
leritil  there  mtiy  he  slight  heave  and  cracking  at  the 
surface.  In  some  mate,  htls  (e.g.  soils  or  weak 
rocks),  the  cavity  can  collapse  and  cause  minor 
subsidence  at  the  ground  surface,  even  when  the 
absolute  values  of  charge  depth  and  charge  si/e 
tire  small.  Most  rocks  collapse  to  form  these  verti¬ 
cal  "chimneys"  when  nuclear  explosions  occur  at 
great  depth. 


Common  Shallow  Croier 


'  Con  tie  SmoU  Lump 
or  Pit  o?  8ose 


Deep  Conicoi  Oo^ef 


Deep  Pounded  Croter 


}  rjurc  21.  (  ross  sci  linns  of  afyparcni  crauns. 


(raters  are  usually  thought  of  as  conical  or  bow  1- 
shaped,  and  as  a  general  idea  this  is  not  far  from 
the  truth,  for  the  true  crater  and  often  for  the  ap¬ 
parent  crater.  However,  apparent  craters  can  vary 
greatly,  from  pits  to  surface  mounds,  depending 
on  charge  depth,  variations  in  ground  material, 
and  differences  in  the  behavior  of  the  fallback  de¬ 
bris.  Figure  21  gives  some  examples  of  cross-section 
shapes  for  apparent  craters. 

Soviet  writers  use  a  ‘‘crater  shape  characteristic” 
n  that  is  the  ratio  of  crater  radius,  r,  to  charge 
depth,  f4i  i  t-'-  ''  =  (‘'Ce  Henrych  1979).  The 

value  n  =  1  (i.e.  r  =  d^)  gives  the  so-called  ‘‘stan¬ 
dard  crater."  As  can  be  seen  from  the  data  com¬ 
pilations  that  are  given  later  in  this  section,  n  =  1 
is  hardly  "standard"  in  rocks  and  other  compe¬ 
tent  materials.  If  this  value  is  reached  at  all,  it  oc¬ 
curs  with  charges  that  are  just  a  bit  shallower  than 
critical  depth. 


CRATKR  .SCALING 

Crater  dimensions,  and  also  the  trajectories  of 
tlyrock,  are  scaled  in  order  to  remove  the  effect  of 
charge  size,  thereby  facilitating  comparisons  of 
test  data.  As  for  other  explosion  phenomena  that 
involve  concentrated  charges  and  three-dimen¬ 
sional  effects,  an  obvious  first  approximation  for 
adjusting  crater  dimensions  is  cube  root  scaling, 
which  implies  geometric  similarity.  In  other 
words,  linear  dimensions  are  scaled  with  respect  to 
charge  radius  (giving  a  dimensionless  variable),  or 
with  respect  to  the  cube  root  of  charge  weight  or 
energy  (giving  physical  dimensions  that  are  ugly 
but  convenient).  However,  both  theoretical  rea¬ 
soning  and  experimental  results  suggest  that  there 
are  other  possibilities  f  jr  scaling  relations. 

To  simplify  the  considerations,  we  first  assume 
that  the  ground  material  is  isotropic,  and  that 
there  is  no  overall  variation  of  bulk  properties 
(such  as  density,  strength,  and  modulus)  with 
depth  or  with  radius  from  the  shot  point.  We  also 
ignore  variation  of  explosive  type,  assuming  that 
the  mass  and  energy  of  the  charge  are  proportion¬ 
al  with  a  constant  ratio. 

If  crater  size  is  controlled  by  the  inherent  strength 
and  density  of  the  material,  the  effects  of  constant 
gravity  can  be  neglected.  Dimensional  analysis 
then  shows  crater  volume  to  be  proportional  to 
either  charge  mass  or  charge  energy  (see,  for  ex¬ 
ample,  Holsapple  and  Schmidt  1980).  For  geomet¬ 
rically  similar  craters,  corresponding  linear  dimen¬ 
sions  are  proportional  to  the  cube  root  of  charge 
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mass,  or  (he  cube  root  of  charge  energy.  LKe  of 
cube  root  scaling  in  rock  blasting  goes  back  to  the 
nineteenth  century,  when  presumably  it  was  de¬ 
duced  empirically. 

By  contrast,  if  gravitational  effects  are  consid¬ 
ered  significant,  the  component  of  material 
strength  that  derives  from  gravity  body  forces  in¬ 
creases  with  depth,  and  a  significant  amount  of 
energy  is  needed  to  lift  debris  from  the  crater 
against  gravity.  If  gravitational  forces  become 
dominant  in  controlling  the  crater  formation,  and 
if  the  gravitational  acceleration  remains  constant, 
dimensional  analysis  shows  crater  volume  to  be 
proportional  to  charge  energy  raised  to  the  power 
'4.  Thus,  for  geometrically  similar  craters,  corre¬ 
sponding  linear  dimensions  would  be  proportional 
to  energy,  or  charge  mass,  raised  to  the  power  'A. 
This  is  called  fourth  root  scaling,  or  gravity  scal¬ 
ing  (White  1973,  Holsapple  and  Schmidt  1980). 
Fourth  root  scaling  also  applies  to  underwater  ex¬ 
plosions  when  gravitational  effects  dominate  over 
hydrodynamic  effects  (see  Part  II).  Referring  to 
Part  II,  this  can  be  appreciated  from  eq  8,  eq  18 
and  Figure  20  when  conditions  are  such  that  the 
charge  deinh  is  proportional  to  the  maximum  bub¬ 
ble  radius  and  is  also  much  greater  than  the  atmos¬ 
pheric  head.  C  omparison  with  underwater  explo¬ 
sions  brings  out  the  point  that,  for  cratering  in  im¬ 
permeable  solids,  the  effective  static  pressure  is 
the  overburden  pressure  /;/t/s  the  atmospheric 
pressure. 

In  general,  it  would  be  reasonable  to  expect 
cube  root  scaling  to  apply  where;  1 )  the  scale  of  the 
problem  is  small,  so  that  depths  in  the  ground  (or 
gravity  body  forces)  are  also  small;  and  2)  the 
ground  material  has  high  inherent  cohesive 
strength  (not  depending  much  on  the  overburden 
pressure).  Fourth-root,  or  gravity,  scaling  might 
look  like  a  better  bet  where:  1 )  very  big  charges  are 
being  used,  so  that  depths  and  overburden  pres¬ 
sures  are  relatively  large;  and  2)  the  ground  mater¬ 
ial  has  little  inherent  cohesion,  deriving  its 
strength  mainly  from  overburden  pressure  (like 
dry  sand). 

Most  materials  have  finite  shear  strength,  even 
under  zero  gravity,  and  in  all  terrestrial  situations 
there  are  finite  gravity  body  forces.  Thus,  in  the 
real  world,  it  might  be  expected  that  the  scaling  ex¬ 
ponent  for  linear  dimensions  would  vary  between 
'A  and  ‘A.  However,  there  are  other  complica¬ 
tions,  both  theoretical  and  practical.  Among  the 
latter  are:  I)  variation  of  effective  rock  strength 
with  the  stressed  volume  (flaw  probability),  2)  var¬ 
iation  of  rock  strength  with  strain  rate  (see  Ciaff- 


ney  1984a, b),  and  3)  the  inevitability  ot  rock  va¬ 
riation  with  depth. 

Direct  experimental  determination  of  crater 
scaling  is  quite  difficult.  Charge  weight,  or  energy, 
has  to  vary  by  several  orders  of  magnitude  to  get 
reliable  results.  At  the  same  time,  scaled  charge 
depth  has  to  be  maintained  constant,  even  though 
the  scaling  relation  for  charge  depth  is  not  known 
with  certainty.  Ideally,  the  ground  material  should 
be  completely  uniform,  so  that  all  tests  are  made 
in  the  same  material.  Finally,  the  same  explosive 
should  be  used  for  all  tests  (equal  energy  for  HE 
and  nuclear  explosions  does  not  necessarily  pro¬ 
duce  the  same  mechanical  effects).  The  experi¬ 
mental  scatter  in  field  data  is  considerable,  and 
replicate  tests  are  needed  to  establish  representa¬ 
tive  values  for  the  crater  dimensions.  However,  it 
is  simply  not  feasible  to  have  replications  of  very 
large  experimental  cratering  explosions,  especially 
nuclear  ones. 

Because  of  these  complications,  test  results  do 
not  always  define  clear  relations  when  crater  di¬ 
mensions  are  plotted  against  charge  weight  on  log¬ 
arithmic  scales.  Some  specialists  take  limited 
available  data  for  particular  ground  materials  and 
particular  moisture  conditions  and  then  derive 
scaling  exponents  by  power  law  regression.  This 
procedure  results  in  exponents  between,  or  even 
outside,  the  limits  of  'A  and  '^4,  and  exponents  are 
often  given  to  three  significant  figures.  Such  a 
procedure  does  not  commend  itself  to  this  writer; 
there  is  no  strong  reason  to  doubt  that  cube  root 
scaling  is  applicable  to  most  small  explosions. 


l  i^iirc  22.  Radius  of  the  apparent  crater  piittted  asiaiust 
ehartte  weight.  Ihtta  are  for  T\ T  e.xplosions  m  dr\  ela  v 
at  a  sealed  de/>th  oj  0.51  It  Ih  tO.2  lu  A  c  j.  tl  roui 
data  eompded  hv  Ronke  et  a!.  1974.) 


Equivolent  TNT  Charge  Weight  (lb) 

/  laurc  2.^.  Ru(/nis  o  f  ific  apparcnf  i.  rutcr  pfoUcd  against  charge  weight  for  chetuical 
and  nuclear  cxplo'.iitns  in  vnrioii'.  f;n>iind  nuiicrials.  (I  rom  data  compiled  hr 
Ro<ikc  Cl  at.  0^4.) 


10'  lo'  (0“  10"  10® 


Equivalent  TNT  Chorge  Weight  (lb) 

Huinv  24.  Radius  of  true  crater  plotted  aiiainst  charge 
wistht  for  r,\  7  c.xplosions  in  sandstone  at  a  sealed  charge 
depth  oj  0.37  ft  Ih  10.15  in  ki;  ).  (I  roin  data  compiled 
hr  Rooke  et  at.  1974.) 


where  gravity  body  forces  are  small,  while  a  slow 
transition  to  fourth  root  scaling  is  to  be  expected 
as  explosion  size  increases.  Figures  22-24  demon¬ 
strate  the  applicability  of  cube  root  scaling  in  dry 
clay  and  hard  rock  for  the  apparent  crater  radius 
from  large  chemical  explosions.  However,  for 
charges  in  the  range  1  to  100  tons,  a  scaling  expo¬ 
nent  of  0.3  seems  to  give  better  representation  of 
the  data  (Fig.  24). 

CRATKR  DIMKNSIONS 

Crater  dimensions  for  a  given  material  are  usu¬ 
ally  presented  as  plots  of  radius,  depth  and  vol¬ 
ume  against  charge  depth.  Linear  dimensions  may 
be  scaled  with  respect  to  the  cube  root  of  charge 
weight,  or  with  respect  to  equivalent  charge  ra¬ 
dius.  Volume  is  usually  scaled  with  respect  to 
charge  weight,  or  energy  yield.  Because  of  the  lim¬ 
itations  of  cube  root  scaling  that  were  discussed 
above,  the  scaling  factor  may  be  charge  weight 
raised  to  an  arbitrary  power  (approximately  0.3). 
Most  of  the  available  data  are  for  the  apparent 
crater,  which  is  easy  to  measure,  although  for 
some  engineering  purposes  the  dimensions  of  the 


true  crater  are  more  informative  and  useful.  When 
the  results  of  many  cratering  shots  with  widely 
varying  charge  weights  and  charge  depths  are 
compiled  for  a  single  ground  material,  there  is 
usually  wide  scatter.  The  scatter  tends  to  be  great¬ 
est  lor  charges  that  are  close  to  optimum  depth, 
probably  because  of  the  abrupt  transition  from 
optimum  depth  to  critical  depth.  In  this  section, 
some  lest  results  are  shown  as  data  bands  which 
represent  the  scatter  that  can  be  expected  due  to 
variations  in  ground  conditions  and  explosive 
type.  Other  tests  results  are  necessarily  shown  as 
simple  curves,  since  that  form  of  presentation  was 
used  by  the  original  investiBators.  However,  in 
reality  test  results  do  not  del.ae  clean  curves  un¬ 
less:  I )  the  test  site  has  highly  uniform  ground  ma¬ 
terial  (no  variation  of  properties,  hori/ontally  or 
vertically),  2)  a  single  explosive  type  is  used,  and  3) 
the  geometry  ol  charge  and  shothole  is  si?nilar  for 
all  shots. 

l  or  HI'  charges  somewhere  around  I  ton.  .lohn- 
son  (1971)  suggested  scaling  with  respect  to  charge 
weight  raised  to  the  power  0.3.  l  igure  25  gives 
scaled  dimensions  for  the  apparent  crater  based  on 
a  I -ton  charge  ol  I'N  1  and  a  scaling  exponent  of 
0.3. 


- D r  ^  Rock 

- nrv  Soti  . Soiuroted  Cloy  Shole 


Ci^orge  Depth  (ft/ton''") 

Figure  25.  Scaled  dimensions  of  apparent 
craters  in  three  kinds  of  ground  materials. 
Scaling  exponent  taken  as  0.3.  (.After  John¬ 
son  1971.) 


Rooke  et  al.  (1974)  used  test  data  and  empirical 
scaling  to  derive  representative  crater  dimensions 
for  a  1-ton  sphere  of  TNT  in  a  wide  variety  of  ma¬ 
terials.  The  results  for  apparent  craters  are  sum¬ 
marized  in  Figure  26.  To  use  these  curves  for  pre¬ 
diction  or  for  blast  design,  it  would  seem  reasona¬ 
ble  to  use  cube  root  scaling  when  scaling  down 
from  1  ton  to  significantly  smaller  charge  sizes, 
and  perhaps  something  closer  to  fourth  root  scal¬ 
ing  when  scaling  up  to  significantly  bigger  HE 
charges. 

In  Figure  27,  dimensions  of  apparent  craters 
made  by  256-lb  (1 16-kg)  charges  of  TNT  in  dry  al¬ 
luvium  are  shown  in  unsealed  form.  The  data 
bands  give  some  indication  of  the  uncertainty  that 
exists  for  a  single  material  with  identical  charges 
of  the  same  explosive.  Cube  root  scaling  should 
give  an  acceptable  approximation  if  these  results 
are  used  to  predict  for  charges  in  the  range  50  to 
1200  lb  (25  to  550  kg). 

A  set  of  data  for  apparent  craters  in  “sandy 
overburden”  blasted  by  ANFO  or  sensitized  slurry 
is  shown  in  Figure  28.  The  charge  weights  were  25 
lb  and  125  lb. 


Figure  26.  Dimensions  of  apparent  craters  formed  in  a  range  of  ground  materials  hy  the  explosion  of  a  I -ton 
charge  of  TNT.  (After  Rooke  et  a!.  1974.) 
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I  mure  27.  Diniension',  oj  the  apparent  erali'r  in  dry  aHu-  l  iipire  2f<.  Sealed  dimensions  of  apparent 

viion  alter  e.xp/osion  of  a  256-lh  H  16-ky)  e/iarye  of  7.\'T.  eralers  in  sandy  overburden.  (After  Bauer  el 

I.  Xfler  Cooper  197b.)  at.  1972.) 


Figure  29.  Dimensions  of  apparent  craters  formed  by  e.xtremely 
small  charges  in  Ottawa  sand.  (After  Piekutowski  1974.) 


Figure  29  gives  cratering  curves  for  blasts  in  dry 
Ottawa  sand  by  tiny,  pea-size  charges  of  lead 
azide.  A  remarkable  thing  about  these  results  is 
that  they  scale  so  as  to  give  credible  values  for 
much  bigger  charges  (Table  3). 

Some  data  on  true  crater  dimensions  for  small 


charges  in  a  variety  of  rocks  are  given  in  Figures 
30-33.  The  results  for  granite  do  not  conform  to 
typical  crater  characteristics,  as  can  be  seen  from 
the  plot  of  crater  depth.  The  charges  in  granite  did 
not  break  out  the  rock  immediately  surrounding 
the  charge  itself;  see  Duvall  and  Atchison  (1957) 
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l  ahle  3.  Approximate  dimensions  of  apparent  eralers  in  rocks  and  soils  (conventional  explosives). 
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for  details.  The  same  type  of  behavior  was  found  Nuclear  explosions  differ  from  chemical  ones — 

for  very  small  charges  in  granite  (D’Andrea  et  al.  energy  density,  temperature  and  pressure  at  the 

1970),  but  the  crater  depth  data  did  not  scale  well.  source  are  all  higher,  and  gas  expansion  is  pro- 

Figure  34  gives  a  general  idea  of  crater  radius  and  duced  indirectly  from  the  surrounding  material  by 

crater  volume  for  very  small  charges  (s  I  lb)  in  vaporization.  This  suggests  that  dimensions  for 

granite.  On  a  larger  scale,  dimensions  of  true  nuclear  craters  might  be  different  than  those  for 

craters  produced  by  1 -ton  charges  in  .sandstone  are  HE  craters,  but  some  of  the  test  data  show  close 

indicated  in  Figure  35.  similarity  between  the  two  types  when  an  appro- 
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f  iMiire  30.  Scaled  dimensions  of  true  craters  in  chalk. 
(After  Duvall  and  Atchison  1957.) 
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Figure  31.  Scaled  dimensions  of  true  craters  in  gran¬ 
ite.  (After  Duvall  and  Atchison  1957.) 


priate  energy  factor  is  applied.  From  considera¬ 
tion  of  theoretical  simulations  and  empirical  com¬ 
parisons,  it  has  been  accepted  by  some  authorities 
that  1  kiloton  of  nuclear  energy  gives  about  the 
same  cratering  performance  as  10‘  lb  of  TNT,  i.e. 
for  crater  comparisons  the  “weight”  of  a  nuclear 
device  is  multiplied  by  0.5.  However,  it  is  not  clear 
to  this  writer  whether  the  factor  is  based  on  direct 
comparisons  for  equal  energy,  or  whether  the  fac¬ 
tor  simply  produces  compatibility  according  to 
some  preconceived  scaling  rule. 
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I  inure  J2.  Sealed  dimensions  oj  true  craters  in  sandstone. 
I  After  Duvall  and  Atchison  1957.) 
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Hfnure  33.  Sealed  dimen  sion  s  of  true  craters  in  marl- 
stone.  (A  fter  Duval!  and  Atchison  1957.) 


A  set  of  curves  for  nuclear  apparent  craters  in 
soil  and  rock  is  given  in  Figure  36.  The  depth 
scales  on  these  curves  imply  applicability  of  cube 
root  scaling.  Representative  dimensions  for  nu¬ 
clear  apparent  craters  in  soils  and  rocks  are  given 
j  in  scaled  form  in  Figure  37. 

Figure  38  gives  dimensions  of  true  craters  from 
j  l-kt  nuclear  bursts  in  rock  and  desert  alluvium. 

For  rock,  cube  root  scaling  is  apparently  a  good 

I 


approximation.  For  alluvium,  cube  root  scaling 
applies  for  crater  depth  and  charge  depth,  but  the 
empirical  scaling  for  radius  is  apparently  closer  to 
fourth  root.  Figure  39  gives  another  curve  for  the 
true  crater  radius  in  hard,  dry  rock.  The  dashed 
line  shows  the  relation  for  true  crater  depth  sug¬ 
gested  by  the  original  author  (Hughes  1968);  this 
type  of  relation  can  apply  only  if  crater  depth  is 
redefined  as  the  depth  below  surface  of  the  base  of 
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l^iire  34.  Scaled  dimensions  of  crater 
onncd  by  small  charges  in  granite. 
A  fter  D'. Andrea  et  al.  1970.) 
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Figure  36.  Dimensions  of  apparent  craters  formed  in 
soils  and  rock  by  l-kiloton  nuclear  explosions.  The 
same  data  are  given  in  scaled. form  by  Hughes  (1968), 
using  a  scaling  exponent  of  i/3.4,  i.e.  0.294.  (After 
Rooke  et  al.  1974.) 
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Figure  35.  Dimensions  of  true  craters 
produced  in  sandstone  by  I -ton  charg¬ 
es  of  T\T.  (After  Rooke  et  at.  1974.) 


the  explosion  cavity,  irrespective  of  whether  it  is  a 
crater  or  a  camouflet  chamber. 

The  cratering  curves  given  in  this  section  and 
the  next  one  represent  test  results  from  large  num¬ 
bers  of  explosions  in  a  wide  range  of  ground  ma¬ 
terials.  At  first  sight  the  curves  may  appear  con¬ 
fusing  or  even  contradictory,  but  when  the  data 
are  compared  in  scaled  form  they  are  very  consis¬ 
tent.  Tables  3  and  4  give  some  characteristic  di¬ 
mensions  for  apparent  craters  and  true  craters  in 
soils  and  rocks  when  the  charge  size  is  up  to  1  ton 
or  so.  Tables  5  and  6  give  comparable  data  for 
craters  formed  by  nuclear  explosions.  It  turns  out 
that  there  is  surprisingly  little  difference  in  the 
crater  dimensions  for  materials  that  have  widely 
different  properties.  For  example,  the  maximum 
scaled  radius  of  the  crater  varies  only  by  a  factor 
of  2  over  the  whole  range  of  soils  and  rocks,  in¬ 
cluding  frozen  soils  and  ice  (see  Tables  10  and  1 1 ). 
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I  3^.  Sailed  (liinensioiis  of  apparent  craters  fanned  in  soils  and  rocks  by  nuclear  explosions. 
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iiyure  38.  Dimensions  of  true  craters  formed  in  rock  and  alluvium  by  J- 
kiloton  nuclear  explosions.  (.After  Rooke  et  a! .  1974.) 


In  order  to  develop  a  practical  feel  for  cratering, 
it  is  useful  to  have  some  easily  appreciated  num¬ 
bers  for  describing  crater  proportions.  In  Tables  7 
and  8,  the  ratio  of  crater  depth  to  crater  radius  is 
listed  for  apparent  craters  and  true  craters,  and  an 
equivalent  sideslope  angle  is  given.  The  numbers 
are  derived  from  cratering  curves  given  in  Figures 
22-24,  and  from  other  sources. 

For  apparent  craters,  the  equivalent  sideslope 
for  most  materials  tends  to  be  around  20"  when 


the  charge  depth  is  zero.  When  the  charge  is  deep 
enough  to  produce  maximum  crater  radius,  the 
equivalent  sideslope  is  around  26°  to  30°.  When 
the  charge  depth  is  such  that  it  gives  maximum 
crater  depth,  the  equivalent  sideslope  is  just  a  bit 
steeper,  say  about  28°  to  32°. 

There  is  very  little  systematic  information  on 
true  craters  in  rock.  For  zero  charge  depth,  the 
first  three  entries  in  Table  8  suggest  20°  to  26°  for 
(he  equivalent  sideslope,  while  the  last  four  entries 
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Table  5.  Approximate  dimensions  of  nuclear  apparent  craters  in  soil 

and  rock. 
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Table  6.  Approximate  dimensions  of  nuclear  true  craters  in  soil  and 

\raximum  crater  radius,  Maxunum  crater  depth, 
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(representing  small-scale  tests)  show  shallower 
slopes.  When  the  charge  is  deep  enough  to  pro¬ 
duce  maximum  radius,  the  equivalent  sideslope  of 
the  true  crater  is  steep  in  most  cases — 36"  to  47  '. 
I  wo  entries  in  Table  8  (small-scale  tests  in  sand¬ 
stone  and  granite)  give  very  shallow  angles  for 
deep  charges,  but  tiiese  results  arc  suspect. 

Although  the  cube  root  scaling  of  linear  crater 
dimensions  is  not  accepted  universally,  there  is 
genera!  agreement  that  crater  volume  can  be 
sealed  with  respect  to  the  charge  weight,  or  the 


energy  yield.  This  is  equivalent  to  assuming  con¬ 
stant  specific  energy  for  the  cratering  process,  irre¬ 
spective  of  the  size  of  the  event.  In  fact,  the  recip¬ 
rocal  of  volume  per  unit  weight  is  a  specific  energy 
if  multiplied  by  the  energy  per  unit  weight  of  the 
explosive.  The  scaled  volume  itself  is  the  recipro¬ 
cal  of  the  traditional  “powder  factor”  used  in 
blasting  practice. 

As  the  charge  depth  increases  from  zero,  crater 
volume  increases  up  to  some  maximum  value, 
after  which  it  falls  again,  reaching  zero  at  the  criti- 
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f  i^urc  39.  Scaled  dimensions  of  true  craters 
formed  in  hard  rock  by  nuclear  e.xplosions. 
(After  Huyhes  1968.) 
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f  igure  40.  “Cratering  efficiency" for  high  e.\- 
plosives  in  various  rocks  and  soils.  (Cooper 
1976.) 
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Figure  41.  Specific  volume  of  the 
apparent  crater  as  a  function  of 
seated  charge  depth  in  soft  rock. 
(Reinterpretation  of  data  from 
Fig.  40.) 


Figure  42.  Specific  volume  of  the 
I  apparent  crater  as  a  function  of 
sealed  charge  depth  in  hard  rock. 
(Reinterpretation  of  data  front 
Fig.  40.) 
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Figure  43.  Specific  volume  of  the  apparent  crater  as  a  function  of  scaled  charge  depth  in  wet 
soils.  (Reinterpretation  of  data  from  Fig.  40.) 


cal  charge  depth.  Figures  27-29  give  examples  of 
the  trend  for  apparent  crater  volume;  Figures  30- 
34  show  the  trend  for  true  crater  volume. 

For  consideration  of  weapons  effects,  the  blast¬ 
er’s  terms  ‘‘powder  factor”  (Ib/yd’)  or  ‘‘specific 
charge”  (kg/m’)  have  been  replaced  by  the  inap¬ 
propriate  term  ‘‘cratering  efficiency,”  typically 
expressed  in  units  of  ft’/ton  for  HE  craters  (which 
is  really  specific  volume).  Since  crater  volume  is  of 
secondary  concern  in  the  context  of  weapons  ef¬ 
fects,  test  data  are  less  readily  available  than  are 
data  for  the  linear  dimensions  of  craters. 

High  explosive  test  data  for  a  variety  of  ground 
materials  have  been  summarized  in  'one  of  the  lit¬ 
erature  by  plotting  the  ‘‘cratering  effiency” 
against  a  dimensionless  charge  depth  (Fig.  40). 
This  dimensionless  charge  depth  is  the  actual 
charge  depth  normalized  with  respect  to  the  cube 
root  of  the  apparent  crater  volume  for  that  charge 
depth.  Such  plots  are  not  directly  useful  for  pre¬ 
diction  and  design,  so  Figures  41  -44  give  the  same 


information  converted  to  conventional  cube  root 
scaling.  The  range  of  the  data  is  insufficient  to 
show  the  decrease  in  specific  volume  as  charge 
depth  approaches  critical  depth,  but  the  results  for 
negative  charge  depth  (air  bursts)  illustrate  how 
crater  volume  becomes  insignificant  once  the 
charge  ceases  to  be  in  contact  with  the  ground  sur¬ 
face.  A  spherical  charge  of  typical  high  explosive 
ceases  to  contact  the  surface  when  its  center  is  at  a 
scaled  height  of  more  than  0.15  ft/lb''''  (0.06  m/ 
kg'''’)- 

When  c/atering  data  for  high  explosives  are 
summarized,  the  standard  reference  explosive  is 
TNT.  The  cratering  performance  for  a  different 
explosive  can  be  accounted  for  by  multiplying  the 
actual  weight  of  that  explosive  by  a  conversion 
factor  in  order  to  obtain  an  approximately  equiva¬ 
lent  weight  of  TNT  (Table  9).  These  factors,  espe¬ 
cially  the  one  for  nuclear  explosions,  should  be  ap¬ 
plied  with  caution. 
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Figure  44.  Specific  volume  of  the  apparent  crater  as  a  function  of  scaled  charge 
depth  in  dry  sod.  (Reinterpretation  of  data  from  Fig.  40.) 


Table  9.  Adjustment  factors  for  cra¬ 
tering  efficiency.  (Non-nuclear  data 
from  Rooke  et  al.  1974.) 

Multiply  weight  of  explosive  by  F  to  obtain 
equivalent  TNT  weight. 


Explosive 

F.  Cratering 
factor 

TNT* 

1.00 

Dynamite  (40<f'o) 

0.68 

Ammonium  nitrate 

1.00 

Nitromethane 

1.10 

C-3,  C-4 

1.34 

Pentolite  (PETN/TNT) 

1.23 

Amatol  (AN/TNT) 

0.94 

Nuclear^ 

0.5 

•  Heat  of  detonation  taken  as  10'  cal/ion  or 
10  cal/kiloton. 

t  Nominal  kiloton  or  megaton  "weight”  to 
be  multiplied  by  F. 


CRATERS  IN  FROZEN  GROUND, 

ICE  AND  SNOW 

Apparent  craters  in  frozen  materials 

Test  data  for  apparent  craters  in  frozen  soils, 
massive  ice,  and  dense  snow  have  been  compiled 
and  plotted  in  appropriate  combinations  (Mellor 
1985).  Representative  data  bands  for  these  plots 
are  shown  in  Figures  45-48,  permitting  easy  com¬ 
parison  with  corresponding  data  for  common 
rocks  and  soils.  Characteristic  crater  dimensions 
are  summarized  in  Table  10. 

Figure  45a  represents  the  combined  data  for 
variation  of  crater  radius  with  charge  depth  in  fro¬ 
zen  silt.  The  general  trend  is  fairly  clear,  with  radi¬ 
us  reaching  its  maximum  value  of  1.8  to  2.7  ft/lb 
(0.7  to  1.1  m/kg'''’)  at  a  scaled  charge  depth 
around  1 .7  to  2. 1  ft/lb (0.7  to  0.8  m/kg  Criti¬ 
cal  charge  depth,  at  which  crater  radius  becomes 
zero,  appears  to  be  around  2.7  to  3.3  ft/lb'''  (0.8  to 
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Figure  46.  Scaled  dimensions  of  apparent  craters  in  frozen  till.  (Data  from  Livingston  and  Murphy  1959,  Mellor  and 
Sellmann  1970,  Bauer  et  al.  1974,  Smith  and  Mellor  1975,  Smith  1980,  Simpson  1981.) 
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Figure  47.  Scaled  dimensions  of  apparent  craters  in  >na.s.sive  ice.  (Data  from  Livingston  I960.) 
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Table  10.  Dimensions  of  apparent  craters  in  frozen  ground,  ice  and  snow. 


Crater  radius,  ft  Ih  ' 

Crater  depth,  fr/lh  ' 

Critical  charge  depth,  ft /lb  ' 

(m/kg 

1 

(m/kg  ■) 

(m'/kg/fur 

Mu.ximuff} 

Optimum 

.Ma.\imum 

Optimum 

Zero 

Zero 

crater 

charge 

crater 

charge 

crater 

crater 

Mufcriu/ 

radius 

depth 

depth 

depth 

radius 

depth 

1  ro/cn  mIc 

1.8-:.? 

1. 7-2.1 

0.6-1. 3 

=  I 

2. 7-3. 3 

2. 7-3. 3 

(0.7-1. 1) 

(0.7-0.8) 

(0.2-0.5) 

(  =  0.4) 

(0,8- 1.3) 

(0.8-1. 3) 

I  ro/en  lill 

1,7-2. 2 

1,3- 1.5 

0.7-1. 1 

1.3- 1.4 

1,8-? 

1. 6-2.0 

(0.7-0.9) 

(0.5-0.6) 

(0.3-0.4) 

(0.5-0.6) 

(0.7-?) 

(0.6-0.8) 

Massive  lee 

2.0-2. 8 

1. 1-1.4 

0.8- 1.4 

0.7-1. 5 

(3.4-4.2) 

(2. 5-4.0) 

(0.8-1. 1) 

(0,4-0.6) 

(0.3-0.6) 

(0.3-0. 6) 

(1.3-1, 7) 

(1.0- 1.6) 

Dense  snow 

2.4-4. 3 

2. 5-3.5 

0.4- 1.5 

1. 0-2.0 

— 

— 

(1.0-1, 7) 

(1.0- 1.4) 

(0.2-0.6) 

(0.4-0.8) 

- 

- 

1.3  m/kg  ').  These  values  are  consistent  with  the 
corresponding  values  for  the  true  crater  (see  next 
section),  but  the  maximum  radius  of  the  true 
crater  is  roughly  40%  greater  than  the  maximum 
radius  of  the  apparent  crater.  The  data  do  not  pro¬ 
vide  values  for  the  radius  of  the  apparent  crater 
that  is  produced  by  a  contact  burst  (charge  depth 
zero). 

Figure  45b  represents  the  combined  data  for  the 
depth  of  the  apparent  crater  in  frozen  silt.  Over 
the  range  of  the  data,  there  is  not  much  systematic 
variation  until  the  critical  charge  depth  is  reached. 
From  the  general  downward  trend  of  the  data 
band,  the  apparent  crater  seems  deepest  at  the  rel¬ 
atively  shallow  charge  depth  of  1  ft/lb (0.4 
m/kg"’). 

Figure  46  represents  the  combined  data  for  fro¬ 
zen  till.  The  range  of  values  for  scaled  charge 
depth  is  too  limited  for  definition  of  a  clear  trend. 
As  in  many  blasting  tests,  the  investigators  were 
interested  in  optimum  blast  design,  to  the  exclu¬ 
sion  of  data  extremes.  Both  crater  radius  and 
crater  depth  seem  to  reach  maximum  values  at 
charge  depths  around  1.3  to  1.5  ft/lb'"'’  (0.5  to  0.6 
m/kg '''’).  This  differs  from  the  optimum  range  for 
the  true  crater  (see  below),  which  is  about  1.8  to 
2.1  ft/lb''’  (0.71  to  0.83  m/kg'^'’).  The  maximum 
radius  of  the  true  crater  is  appreciably  bigger  than 
that  of  the  apparent  crater — by  25%  to  65%.  Crit¬ 
ical  depth  for  the  apparent  crater  seems  to  be 
about  1.6  to  2.0  fl/ib  '''  (0.6  to  0.8  m/kg '''’),  which 
is  shallower  than  the  critical  depth  for  the  true 
crater. 

The  original  radius  data  for  apparent  craters  in 
massive  ice  show  very  wide  scatter  as  charge  depth 
approaches  the  critical  value.  Taking  the  data 
bands  for  crater  radius  and  crater  depth  together 


(Fig.  47),  it  looks  as  if  the  apparent  crater  reaches 
maximum  size  when  charge  depth  is  in  the  range 
1.0  to  1.5  ft/lb'^'  (0.4  to  0.6  m/kg  ''’).  By  contrast, 
the  true  crater  (Fig.  51)  reaches  maximum  size 
when  the  charge  depth  is  in  the  range  3  to  4  ft/lb  '^' 
(1.2  to  1.6  m/kg'''’).  For  the  apparent  crater,  criti¬ 
cal  charge  depth  seems  to  be  close  to  3  to  4  ft/lb’  ’ 
(1.2  to  1.6  m/kg'''’). 

The  data  for  other  frozen  materials  are  insuffi¬ 
cient  to  establish  trends.  For  practical  purposes, 
the  results  for  frozen  clay  shale,  and  for  ice-rich 
till  covered  by  peat,  are  not  significantly  different 
from  the  results  for  frozen  silt. 

The  limited  amount  of  data  for  unfrozen 
(thawed)  soil  provides  some  indication  of  the  ef¬ 
fect  of  freezing  on  crater  size.  The  results  for  un¬ 
frozen  silt  show  both  crater  radius  and  crater 
depth  somewhat  greater  than  corresponding  val¬ 
ues  for  the  frozen  state.  The  same  is  true  when  the 
data  for  unfrozen  gravel  are  compared  with  results 
for  frozen  till.  The  actual  data  can  be  found  else¬ 
where;  a  summary  is  given  by  Mellor  (1985). 

Data  for  apparent  craters  in  deep,  dense  snow 
are  shown  in  Figure  48.  Crater  radius  has  maxi¬ 
mum  values  in  the  range  2.4  to  4.3  ft/lb  ''*  (1.0  to 
1.7  m/kg'''),  and  these  values  are  obtained  when 
the  charge  depth  is  in  the  range  2.5  to  3.5  fl/lb'^' 
(1.0  to  1.4  m/kg '''’)•  There  is  not  much  systematic 
variation  of  crater  depth  with  charge  depth  for 
charge  depths  less  than  3  ft/lb (1.2  m/kg' ')■  The 
maximum  depth  of  the  apparent  crater  is  in  the 
range  0.4  to  1.5  ft/lb'''’  (0.2  to  0.6  m/kg'’)  and 
these  maximum  values  occur  when  charge  depth  is 
in  the  range  1  to  2  ft/lb (0.4  to  0.8  m/kg' ').  Ex¬ 
isting  data  for  snow  do  not  provide  a  clear  indica¬ 
tion  of  critical  charge  depth. 
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Figure  48.  Scaled  dimensions  of  apparent  craters  in  dense  .snow.  (Data  from  Fuchs  1957,  Livingston  1968. j 


True  craters  in  frozen  materials 

Measurements  on  true  craters  are  more  difficult 
to  make  than  are  corresponding  measurements  on 
apparent  craters.  Nevertheless,  a  considerable 
body  of  data  exists  for  frozen  silt,  frozen  gravel, 
massive  ice,  and  dense  snow  (Fuchs  1957,  Living¬ 
ston  and  Murphy  1979,  Livingston  1960,  1968, 
McCoy  1965,  Mellor  and  Sellmann  1970,  Mellor 
1971),  These  results  have  been  compiled  and  plot¬ 
ted  elsewhere  (Mellor  1972,  1973).  For  present 
purposes  the  detailed  data  plots  are  replaced  by 
representative  data  bands,  as  shown  in  Figures 
49-52.  Characteristic  crater  dimensions  are  sum¬ 
marized  in  Table  1 1 . 

Figure  49a  gives  an  impression  of  the  variation 
of  crater  radius  with  charge  depth  for  true  craters 


in  frozen  silt.  Radius  reaches  maximum  values 
when  charge  depth  is  in  the  range  2.1  to  2.7  ft/lb  ‘ 
(0.8  to  1.1  m/kg'"').  With  charge  depth  in  this 
range,  the  crater  radius  is  between  2.4  and  4.0  ft/ 
lb"  (1 .0  to  1.6  m/kg"). 

Figure  49b  shows  how  the  depth  of  the  true 
crater  varies  with  charge  depth  in  frozen  silt.  Until 
optimum  charge  depth  is  reached,  the  bottom  of 
the  crater  ranges  from  charge  depth  to  about  0.8 
ft/lb"  (0.3  m/kg")  below  the  charge  depth. 

The  specific  volume  of  the  true  crater  in  frozen 
silt  (Fig.  49c)  indicates  a  well  defined  value  for  op¬ 
timum  charge  depth — about  2.2  ft/lb"  (0.9 
m/kg").  However,  at  this  optimum  charge  depth 
the  range  of  specific  volume  is  wide — from  18  to 
54  ftVlb  (1.1  to  3.4  mVkg). 


Table  11.  Dimensions  of  true  craters  in  frozen  ground,  ice  and  snow. 


Crater  radius,  fl/lh 


Crater  depth,  ft/th  '  Criticalchar^edeplh.ft/lh  ' 


tm/kg  J 


tm'/kg)  for 


Crater  votume, _ 

Maximum  Optimum 


Critical  charge 
depth,  zero 


Maximum 

Optimum 

Ma.ximum 

Optimum 

Zero 

Zero 

crater  sot. 

charge  depth 

specific  vote 

crater 

charge 

crater 

charge 

crater 

crater 

ft'/lbtm'/kg)  fl/lb  tm/kg  ') 

ft /lb 

Uaterial 

radius 

depth 

depth 

depth 

radius 

depth 

volume 

depth 

(m/kg 

f  ro/on  silt 

4-4.0 

2.1  2.7 

2. 0-3. 2 

2. 3-3.0 

2.9-? 

2. 8-3. 2 

18-54 

=  2.2 

2.8-7 

(1.0-1. 5) 

(O.H-1.1) 

(O.S-1.3) 

(0.9- 1.2) 

(1. 2?) 

(1,1-13) 

(1.1 -3. 4) 

(  =0.9) 

(l.l-'’l 

I To/cn  (ill 

:  1  4.: 

l.y-2.7 

!.ft-3.4 

2.0-2. K 

2. 2-3. 2 

2.2-3, 1 

13-35 

2. 1-2. 3 

2.2-2  9 

(O.S  l.^l 

(O.S  1.1) 

(O.ft-1.3) 

lO.K-l.l) 

(0.9- ), 3) 

(0.9- 1.2) 

(0.8-2.21 
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Massive  leu 

VI  -4. y 
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_ 

— 
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_ 
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- 
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~ 
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_ 

— 
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— 
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— 
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Figure  49.  Scaled  dimensions  of  true  craters  in  frozen  silt. 
(Data  from  McCov  1965,  Mellor  and  Sellmann  1970,  Mellor 
1971.) 


Figure  50a  gives  the  general  magnitude  of  true 
crater  radius  in  frozen  till.  The  data  used  for  this 
plot  were  obtained  in  frozen  tills  which  contained 
both  gravel  and  fine-grained  soil,  and  also  in  clean 
ice-bonded  gravels.  The  upper  limit  of  crater  ra¬ 
dius  in  frozen  till  is  not  much  different  from  the 
upper  limit  in  frozen  silt,  and  optimum  charge 


depths  are  similar  for  the  two  materials.  However, 
because  some  tills  are  stronger  and  denser  than 
typical  silts,  the  lower  limit  of  crater  radius  is  low¬ 
er  for  tills  than  for  silt,  and  critical  charge  depth 
can  be  smaller  in  tills  than  in  silt.  Crater  radius  has 
maximum  values  in  the  range  2.1  to  4.2  ft/lb (0.8 
to  1.7  m/kg''''),  and  these  maximum  values  are  ob- 


( ft/lb''’) 
(m/kg''^) 
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Figure  50.  Scaled  dimensions  of  true  craters  in  frozen  till  or 
frozen  gravel.  (Data from  Livingston  and  Murphy  1959,  Mellor 
and  Sellmann  1970,  Smith  and  Mellor  1975.) 


tained  when  charge  depth  is  in  the  range  i  .9  to  2.7 
ft/lb (0.8  to  l.I  m/kg'^'). 

Figure  50b  represents  the  variation  of  crater 
depth  with  charge  depth  in  frozen  till.  The  trend  is 
broadly  similar  to  that  for  frozen  silt,  but  crater 
depth  is  smaller  over  most  of  the  range  of  charge 
depths,  and  critical  depths  are  somewhat  shallow¬ 
er  in  frozen  till. 

Figure  50c  shows  how  specific  volume  varies 
with  charge  depth  for  true  craters  in  frozen  till. 


Maximum  specific  volume  is  in  the  range  13  to  35 
ftVlb  (0.8  to  2.2  mVkg)  when  charge  depth  is 
close  to  optimum,  i.e.  2.1  to  2.3  ft/lb'  '  (0.8  to  0.9 
m/kg'">). 

True  craters  in  massive  ice  are  appreciably  big¬ 
ger  than  those  in  frozen  ground.  Figure  51a  gives 
crater  radius  as  a  function  of  charge  depth  in  ice. 
Maximum  values  of  radius  are  in  the  range  3.1  to 
4.9  ft/lb  '  (1 .2  to  1 .9  m/kg'  ')  when  charge  depth 
is  in  the  range  3.5  to  4  ftlb' '  (1 .4  to  1.6  m/kg  '). 
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Figure  51.  Scaled  dimensions  of  true  craters  in  massive  ice.  (Data  from  Livingston  I960.) 


Figure  51b  shows  true  crater  depth  in  ice.  The 
form  of  presentation  differs  from  that  used  for 
frozen  ground,  in  that  the  base  of  the  broken  ice  is 
plotted  for  charge  depths  that  exceed  critical  depth 
(the  camouflet  range).  The  depth  of  the  true  crater 
ranges  from  the  base  of  the  charge  to  a  level  that  is 
1.5  to  0.5  ft/lb (0.6  to  0.2  m/kg''^’)  below  the 
base  of  the  charge. 

The  specific  volume  for  true  craters  in  ice  (Fig. 
51c)  has  maximum  values  in  the  range  38  to  118 
ftVlb  (2.4  to  7.4  mVkg)  when  charge  depth  is  in 
the  range  3.3  to  4.2  ft/lb'^’  (1.3  to  1.7  m/kg  '"’)- 
Dense  snow,  such  as  is  found  in  the  surface  lay¬ 
ers  of  polar  ice  caps  and  in  well-settled  seasonal 
snow  packs,  is  much  weaker  than  ice,  and  less 


dense  by  a  factor  of  2  to  3.  Nevertheless,  the  di¬ 
mensions  of  the  true  crater  in  snow  are  not  very 
much  different  from  those  in  ice. 

Figure  52a  shows  crater  radius  as  a  function  of 
charge  depth  for  dense  snow.  Compared  with  cor¬ 
responding  plots  for  ice  and  frozen  ground,  the 
transition  from  optimum  depth  to  critical  depth  is 
less  abrupt.  Maximum  values  of  radius  are  in  the 
range  2.8  to  4.4  ft/lb'''’  (1.1  to  1.7  m/kg'''')  when 
charge  depth  is  in  the  range  2.0  to  2.6  ft/lb (0.8 
to  1 .0  m/kg'''’). 

The  depth  of  the  true  crater  in  snow  is  consis¬ 
tently  below  the  base  of  the  charge — in  the  range 
0.35  to  1.4  ft/lb (0.14  to  0.57  m/kg  '"')  below  the 
charge  base. 
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The  specific  volume  in  snow  has  maximum  val¬ 
ues  in  the  range  45  to  120  ft '/lb  (2.8  to  7.5  m  Vkg). 
These  values  are  obtained  when  charge  depth  is  in 
the  range  3.3  to  4.3  ft/lb  '  (1.3  to  1.7  m/kg'^’). 


ROW  CRATERS,  DITCHES 
A.ND  SINKI.NG  CUTS 

In  military  engineering,  a  line  of  overlapping 
crateis  may  be  requireu  in  order  to  form  an  obsta¬ 
cle.  The  composite  crater  is  usually  termed  a  row 
crater.  In  construction,  explosives  may  be  used 
either  to  form  a  ditch  directly,  or  to  fracture  the 
ground  so  that  a  ditch  can  be  excavated  mechanic¬ 
ally.  The  operation  is  usually  referred  to  as  trench 
blasting,  or  ditch  blasting. 

The  simplest  way  to  make  a  row  crater  or  ditch 
is  to  place  a  line  of  closely  spaced  cratering  charg¬ 
es  at  optimum  depth.  If  the  charges  are  widely 
spaced,  each  produces  its  own  independent  crater 
of  radius  Ri,  following  th^  relations  that  were 
given  previously.  When  the  center-to-center  spac¬ 
ing  of  the  charges  (s)  is  such  that  5  >  23?,,  there  is 
no  overlap  and  the  craters  are  independent,  except 
for  possible  interference  of  the  crater  lips  and 
some  exchanges  of  debris.  However,  values  of  5  up 
to  3/?,  to  4R,  may  be  adequate  for  creating  a  vehi¬ 
cle  barrier  in  military  operations,  since  the  crater 
lips  form  mounds.  Much  closer  spacing  is  needed 
to  develop  a  continuous  row  crater,  or  ditch,  but 
when  s  is  sufficiently  small,  say  s  <  1.43?,,  the 
width  of  the  resulting  ditch  can  exceed  23?,  and  the 
depth  of  the  ditch  can  exceed  the  depth  of  a  single 
crater. 

The  width  and  depth  of  a  row  crater  are  often 
discussed  in  terms  of  an  “enhancement  factor.”  If 


3?,  is  the  radius  of  the  apparent  crater  for  a  single 
charge,  and  D,  is  the  corresponding  depth,  the  ap¬ 
parent  width  of  the  row  crater  and  the  apparent 
depth  of  the  row  crater  Dj  are  such  that  B^/2R,  >  I 
and  Df/D,  >  1  when  s/23?,  <  0.7  (i.e.  s  <  1.43?,). 
The  general  trend  given  by  test  data  is  shown  in 
Figure  53;  the  enhancement  factor  can  reach 
about  1.6  when  s/R,  is  decreased  to  about  0.55. 
To  obtain  optimum  enhancement  in  this  way,  the 
charge  depth  should  be  increased  as  s/R,  decreas¬ 
es.  For  a  given  charge  spacing  s/3?,,  the  optimum 
charge  depth  is  obtained  from  Figure  53  as  d^.  = 
(Bf2R,)d^^  =  (DfD,)dc,,  where  r3(-,  the  charge 
depth  optimized  for  a  single  apparent  crater. 

When  trenches  are  to  be  blasted  for  construc¬ 
tion  work,  various  methods  can  be  used.  If  indi¬ 
vidual  explosive  charges  are  intended  to  break  the 
ground  material  and  to  excavate  it  in  a  single  oper¬ 
ation,  the  procedure  is  identical  to  that  described 
above  for  row  craters;  the  linear  apparent  crater  is 
accepted  as  the  finished  trench.  If  the  ditch  is  for 
land  drainage,  water  flow  can  sometimes  flush  out 
blast  debris  that  is  loose  and  fine-grained.  By  con¬ 
trast,  if  the  operation  is  to  be  more  conservative, 
with  limits  on  explosive  consumption,  flyrock 
(ejecta),  ground  motion,  airblast,  and  overbreak, 
then  the  emphasis  is  on  production  of  overlapping 
true  craters,  with  optimized  breakage.  The  design 
parameters  for  the  true  crater  are  obtained  for  a 
single  crater  that  has  its  charge  depth  just  less  than 
critical  depth,  and  the  degree  of  overlap  is  decided 
by  adjustment  of  both  charge  spacing  and  charge 
depth.  The  broad  aim  is  to  break  up  the  rock  with¬ 
in  the  planned  limits  of  the  excavation,  while 
avoiding  overbreak  and  ejection  of  missiles.  Frag¬ 
mented  material  is  dug  out  with  a  backhoe  or 
something  similar. 


Chorge  Spacing  (s/R,) 

f  igure  53.  Enhancement  of  dimensions  for  the  apparent  crater 
when  closely  spaced  charges  are  fired  in  a  row.  (After  Rooke  et  al. 
1974,  U.S.  Army  1984.) 
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A  single  row  of  cratering  charges  is  less  than 
ideal  lot  excavation  of  a  steep-sided  trench  in 
rock,  but  the  alternative  of  using  multiple  rows  of 
smaller  charges  is  sometimes  prohibitive  in  terms 
of  drilling  and  loading  effort.  One  possibility  is  to 
use  a  single  row  of  vertically  delayed  deck  charges 
(Mellor  and  Sellmann  1970). 

For  pipeline  trenching  in  rock,  the  end  of  the 
trench  may  be  advanced  by  a  technique  similar  to 
bench  blasting. 

Special  explosives  and  delay  detonators  have 
been  developed  for  trenching;  both  the  cartridges 
and  the  caps  resist  premature  detonation  by  prop¬ 
agation. 

.Another  way  to  make  a  large  trench,  or  an  anti¬ 
tank  ditch,  is  to  lay  a  linear  charge  (sausage 
charge,  or  pipe  charge)  in  a  narrow  excavated  cut, 
with  some  backfill  (stemming).  In  this  case,  the 
charge  is  emplaced  like  an  underground  cable  or 
small  diameter  pipe,  using  slot-cutting  equipment 
instead  of  shothole  drills.  Rough  estimates  of  the 
required  charge  per  unit  length  can  be  made  from 
row  crater  data;  if  is  the  required  weight  of 
each  charge  in  a  row  with  charge  spacing  x,  the 
equivalent  charge  per  unit  length  is  U’/s.  Such  an 
estimate  is  only  a  first  approximation,  since  shock 
propagation  and  attenuation  are  two-dimensional 
rather  than  three-dimensional,  but  it  permits  the 


use  of  available  data  from  single-crater  tests.  An 
alternative  approach  is  to  regard  the  buried  linear 
charge  as  a  column  charge  of  the  type  that  would 
be  used  in  bench  blasting.  With  burial  depth  Band 
charge  diameter  D,  the  effective  range  of  B/D  in 
hard  material  might  be  20  to  40. 

Ditches  in  soft  ground  can  be  blasted  by  using 
small  cartridges  of  ditching  dynamite  or  a  modern 
AN  water  gel  equivalent.  The  necessary  shotholes 
can  be  drilled  easily  with  a  simple  auger  or  a  driv¬ 
ing  bar  (punch).  When  cartridges  of  approximate¬ 
ly  1.25  in.  (32  mm)  diameter  are  used,  each  hole 
has  a  collar  distance  of  6  to  12  in.  (150  to  300 
mm),  preferably  fully  stemmed.  Tables  12  and  13 
give  guidelines  for  soft  ground  ditching  with  either 
a  single  row  of  charges  or  multiple  rows.  One  of 
the  traditional  methods  for  ditching  in  wet  soils  re¬ 
lies  on  propagation  (sympathetic  detonation)  be¬ 
tween  the  individual  charges.  Cartridges  of  sensi¬ 
tive  explosive,  usually  ditching  dynamite  (i.e.  50®^o 
straight  dynamite),  are  tamped  into  saturated  soil 
at  fairly  close  spacing,  without  any  kind  of  con¬ 
nection  between  the  individual  charges.  The  first 
charge  in  the  row  is  fired  by  a  single  cap,  and  the 
blast  then  propagates  from  charge  to  charge  by 
shock  transmission.  In  wet  ground,  with  tempera¬ 
tures  above  freezing,  reliable  propagation  can  be 
achieved  over  scaled  distances  of  approximately  2 


Table  12,  Ditching  in  soft  ground  with  single  row  of  charges.  (After 
CM, I..  1984.) 


Cartridges 
per  hole 
!  ‘  4  •  S  in.) 

Depth  to 
top  of 
charge 
(in.) 

Distance 

between 

holes 

(in.) 

Probable 
depth  of 
ditch 

<fi) 

Probable 
lop  width 
of  dill  it 

i.ri> 

P.xplosive 
consumption 
(lb  1(H)  ft) 

,  ^ 

6-8 

12 

1.5 

4 

2s 

1 

6-12 

15 

2.5 

6 

40 

2 

6-12 

18 

3 

8 

65 

3 

6-12 

21 

4 

10 

85 

4 

6-12 

24 

5 

1.1 

1(X) 

S 

6-12 

24 

6 

16 

125 

Table  13.  Ditching  in  soft  ground  with  multiple  rows  of  charges.  (After  CM. I,.  1984b.) 


Cartridges 

Hole 

Row 

Depth  of 

per  hole 

spacing 

'ipacing 

ditch 

4  in.  •  h  in.) 

tin.) 

tin.) 

_  (ft) 

1 

15 

10 

2.5-3 

A 

18 

16 

3  3.5 

.1 

21 

42 

4-4.5 

4 

24 

48 

5-5.5 

5 

24 

48 

6  6.5 

Number  of  rows 


J 

5 

7 

(fii 

(lb/ 1(H)  ft) 

Tftr 

db  l(X)fil 

r/0  (lb  loo.K) 

11 

80 

— 

—  — 

11 

111 

14 

2(X) 

— 

11 

172 

17 

257 

20  141 

17 

200 

21 

KXl 

25  4(X) 

20 

250 

24 

115 

28  5(X) 

9 

ifn  (lb 


I! 

HU  lib  lOOfU 
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40U 

MX) 


29 


5IH) 


M) 


Delay  Period  ^ 
(x  25ms) 


Figure  54.  Example  of  the  layout,  connections  and  delays 
for  charges  used  to  make  a  sinking  cut. 

ft/lb  '  (0.8  m/kg' ').  The  loads  and  spacings  given 
in  Tables  12  and  13  should  provide  reliable  propa¬ 
gation  with  ditching  dynamite  in  wet  soil. 

A  sinking  cut  is  a  multi-charge  blast  that  produc¬ 
es  a  pit  in  a  surface  which  is  close  to  horizontal.  In 
essence,  a  sinking  cut  is  an  array  of  many  cratering 
charges  set  deep  enough  to  avoid  much  throw  or 
flyrock.  Ideally,  the  broken  rock  forms  a  mound 
of  fragments  that  can  be  excavated  easily  by  me¬ 
chanical  equipment. 


The  simplest  form  of  sinking  cut  would  be  an 
array  of  compact  cratering  charges  set  just  shal¬ 
lower  than  critical  depth  and  fired  simultaneously 
to  give  overlapping  craters.  However,  in  typical 
practice  extended  column  charges  in  small-diam¬ 
eter  bore'ioles  are  likely  to  be  used,  and  delays  are 
employed,  both  to  improve  the  breakage  and  to 
reduce  ground  vibrations.  This  makes  the  sinking 
cut  more  like  a  tunnel-blasting  round.  Figure  54 
shows  the  general  arrangement  for  a  sinking  cut 
when  the  pit  is  to  be  roughly  equidimensional  in 
plan.  If  the  pit  is  to  be  long  and  narrow,  more  like 
a  ditch,  then  the  delay  pattern  can  be  arranged  to 
displace  material  towards  the  centerline. 

The  powder  factor  for  sinking  cuts  is  typically 
about  0.75  to  1.5  Ib/yd'  (0  44  to  0.89  kg/m')-  In 
terms  of  specific  volume,  these  values  translate  to 
a  typical  range  of  18  to  36  ftVlb  (1.1  to2.2mVkg), 
with  a  lower  limit  around  11  ftVlb  (0.7  mVkg). 
Such  values  are  not  much  different  from  specific 
volumes  for  the  true  crater  when  a  single  charge  is 
fired  at  optimum  depth  in  common  rocks  (see  Fig. 
30-35). 

Sinking  cuts  are  typically  up  to  15  to  20  ft  (4.6 
to  6.1  m)  deep,  but  can  be  as  much  as  40  ft  (12  m). 
Shotholes  are  usually  drilled  below  the  required 
finished  grade  to  allow  for  the  uneven  surface  that 
results  from  cra-er  blasting. 

Repetitive  crater  blasts  produce  some  enlarge¬ 
ment  of  a  crater.  The  idea  is  to  first  fire  a  charge  at 
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Eiiiure  55.  Effect  of  repeated  explosions  on  the  dimensions  of  apparent  craters  in  rock 
and  soil.  After  Rooke  et  at.  1974,  U.S.  Army  1984.) 
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or  near  the  surface,  then  to  fire  a  second  charge  of 
equal  si^e  at  the  bottom  of  thr  crater,  and  so  on. 
This  procedure,  which  is  referred  to  as  nail  driv- 
inf’,  gives  significant  increase  of  crater  depth  while 
producing  relatively  little  increase  of  the  crater  ra¬ 
dius  at  ground  level.  Figure  55  indicates  the  mag¬ 
nitude  of  the  effect  for  the  radius  and  depth  of 
true  craters  in  rock  and  soil. 


UNDERWATER  CRATERS 
AND  TRENCHES 

In  considering  the  process  of  cratering  under¬ 
water,  it  is  helpful  to  draw  a  parallel  with  standard 
crater  blasting  on  dry  land.  In  one  case  the  soil  or 
rock  is  overlain  with  air,  in  the  other  case  with 
water,  which  is  a  thousand  times  more  dense.  Sys¬ 
tematic  study  of  underwater  cratering  has  not 
gone  very  far,  so  it  is  necessary  to  speculate  on 
how  the  various  processes  might  be  affected. 

For  comparison  of  cratering  on  dry  land  and 
underwater,  we  assume  that  the  water  depth  is  suf¬ 
ficient  to  prevent  direct  blowout  through  the 
water/air  surface.  We  also  recognize  that  the 
strength  S  of  soil  or  rock  has  two  components,  an 
inherent  cohesion  c  and  an  “internal  friction” 
that  increases  with  the  applied  compressive  stress 
p,  such  that  S  =  c  -t-  ptan^,  where  <t>  is  an  effec¬ 
tive  angle  of  internal  friction.  We  also  note  that 
some  soils  and  rocks  are  impermeable,  while 
others  are  permeable,  allowing  pore  water  pres¬ 
sure  to  equilibrate  with  the  surface  fluid  pressure. 

With  these  things  in  mind,  the  following  effects 
seem  likely: 


1.  For  impermeable  underwater  materials,  the 
“internal  friction”  component  of  strength  will  in¬ 
crease  with  increasing  water  depth,  since  p  increas¬ 
es,  but  the  effect  is  small. 

2.  For  permeable  granular  soils  which  have  low 
cohesion  (c  — 0),  submersion  will  reduce  the 
strength  below  the  “dry  land”  value,  since  the  ef¬ 
fective  unit  weight  of  the  soil  grains  is  decreased 
by  buoyancy  (p  is  lower  for  a  given  depth  in  the 
soil). 

3.  When  a  crater  is  formed,  fragmented  mater¬ 
ial  will  settle  to  angles  of  repose  smaller  than  cor¬ 
responding  angles  in  air,  and  loose  material  is  like¬ 
ly  to  wash  back  into  the  crater  by  water  motion. 

4.  When  a  charge  detonates  at  'he  interface  be¬ 
tween  the  water  and  the  bed  material,  the  overly¬ 
ing  water  confines  the  explosion  more  than  air 
would  do  for  a  surface  charge  on  dry  land.  In  ef¬ 
fect,  the  charge  is  mud-capped  by  the  water,  and 
energy  transfer  to  the  underlying  solid  is  likely  to 
be  more  efficient  than  it  would  be  in  air. 

5.  The  breakout  and  heave  of  material  from  a 
buried  explosion  will  be  subdued  by  the  inertial  re¬ 
sistance  of  overlying  water. 

6.  Ejected  fragments  meet  much  higher  resis¬ 
tance  in  water  than  in  air,  since  the  fluid  density  is 
about  a  thousand  times  higher. 

Taking  into  account  all  of  these  potential  ef¬ 
fects,  certain  net  results  can  be  expected.  For 
strong  materials  and  relatively  shallow  water,  a 
charge  set  directly  on  the  bed  could  break  more  of 
the  bottom  material  than  an  equivalent  surface 
charge  in  air  would  do,  but  ejected  fragments  will 
not  travel  very  far  or  very  fast.  A  charge  buried  in 
strong  material  is  likely  to  be  less  effective  than  an 


f  ie’ire  56.  Dimensions  oj  underwater  apparent  <raler\  in  fine-eralned  soils 
as  a  function  of  scaled  water  depth.  One-ton  charge  fired  in  direct  contact 
with  the  bed  of  the  water  body.  (Rooke  et  al.  1974./ 
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equivalent  cratering  charge  on  dry  land,  as  bubble 
expansion  effects  are  subdued  by  overlying  water. 
In  weak  bed  sediments  at  relatively  shallow  depths, 
a  charge  set  directly  on  the  bed  is  likely  to  displace 
more  material  than  an  equivalent  charge  in  air, 
but  the  final  dimensions  of  the  crater  will  be  af¬ 
fected  by  “washback”  and  by  the  shallower  angles 
of  repose  that  are  common  underwater.  Buried 
charges  in  weak  underwater  sediments  could  dis¬ 
place  more  material  than  equivalent  charges  on 
land,  since  the  effective  strength  of  the  soil  is  low¬ 
er  underwater.  However,  the  throw  of  the  material 
is  likely  to  be  suppressed,  leaving  a  relatively  small 
apparent  crater. 

Actual  data  are  sparse.  Figure  56  shows  the  ra¬ 
dius  and  depth  of  the  apparent  crater  for  a  1-ton 
charge  set  directly  on  fine-grained  underwater  ma¬ 
terial,  with  scaled  water  depth  as  a  variable.  At 
very  shallow  water  depth,  the  explosion  vents  di¬ 
rectly  to  air  and  the  crater  dimensions  are  affected 
strongly  by  washback  (unless  the  water  depth  is 
less  than  the  minimum  height  of  the  crater  lip). 
When  the  scaled  water  depth  exceeds  0.8  ft/lb'''’, 
the  crater  dimensions  are  more  than  twice  the  di¬ 
mensions  that  would  be  expected  for  fine-grained 
dry  soils  on  land  (radius  0.9  to  1.3  ft/lb'''*,  depth 
0.3  to  0.5  ft/lb'''’).  There  does  not  appear  to  be  any 
published  information  on  point-charge  craters  in 
rock,  or  on  true  craters  (except  for  situations 
where  the  scaled  water  depth  is  very  small,  in 
which  case  there  is  not  much  difference  from 
blasting  in  saturated  soil  or  rock  on  land). 

Guidelines  for  underwater  trenching  in  coral 
have  been  given  by  Hallanger  (1976),  although  the 
information  provided  is  incomplete.  Diagrams  in¬ 
dicate  the  trench  cross  sections  produced  by  vari¬ 
ous  types  and  amounts  of  explosive  (Fig.  57),  and 
it  is  assumed  here  that  these  cross  sections  are  for 
the  apparent  linear  crater.  Water  depth  is  not 
specified,  but  lack  of  geometric  similarity  as 
charge  weight  increases  suggests  that  the  water  is 
shallow  enough  for  venting  and  washback  to  be 
significant  and  variable  with  charge  size.  The 
charge  per  unit  length  is  known  for  Bangalore  tor¬ 
pedoes  (Fig.  57a);  it  is  approximately  2  Ib/ft  for  a 
single  tube  (9  lb  of  main  charge  plus  a  booster  for 
each  5-ft  length).  Figure  57b  refers  to  packs,  each 
containing  20  lb  of  explosive,  but  the  spacing  is 
not  given  explicitly  (spacing  of  3-4  ft  is  men¬ 
tioned,  perhaps  for  each  heap  of  explosive  pack¬ 
ages,  or  perhaps  only  for  single  20-lb  packages). 
Figure  57b  refers  to  a  non-standard  hose  charge 
designated  “M8  demolition  charge”;  the  results 
are  not  consistent  with  those  for  Bangalore  torpe- 


Trench  Width  (tt) 


Figure  57.  Dimensions  of  underwater  trenches  in  coral 
when  linear  charges  are  fired  directly  on  the  bed: 
a)  Bangalore  torpedoes  (2  Ib/ft),  b)  packs  of  explosive 
(20  lb  each)  at  unspecified  intervals  (perhaps  3-4  ft), 
c)  hose  charges  (unit  weight  not  specified).  (Hallanger 
1976.) 


does,  possibly  because  the  hose  charges  covered 
more  horizontal  width. 

When  multiple  charges  are  used  under  water, 
propagation  is  possible.  This  can  be  an  advantage 
in  providing  a  simple  way  to  initiate  multiple 
charges  of  sensitive  explosive.  It  can  also  be  a  dis¬ 
advantage  if  there  is  a  requirement  for  delays  be¬ 
tween  charges.  To  avoid  propagation,  it  may  be 
necessary  to  use  insensitive  explosives  and  caps. 


BOMB  AND  SHELL  CRATERS 

Representative  dimensions  for  craters  formed  in 
soils  by  old-style  bombs  and  shells  are  provided  by 
measurements  made  during  World  War  II.  Table 
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Table  14.  Representative  dimensions  of  bomb  craters.  (After U.S. 
Army  1965.) 


Oomh 

fkftt 

Depth  of 
burial 

do 

Crater 

diameter 

<fo 

Crater 

depth 

tfo 

Soil 

Remarks 

50 

1.0 

16.5 

4.7 

Clay 

7.1 

19,. 7 

5.2 

Clay 

11.4 

6.0 

2.5 

Clialk 

'» 

8.5 

2.0 

Railway  embankment 

16.0 

9.5 

1.7 

Clay  j 

Deep 

•> 

10.0 

1.0 

Sand  and  gravel 

camouflet 

ItX) 

7.5 

18.5 

5.2 

Clay 

250 

8.0 

53.8 

9.0 

Clay 

•> 

12.0 

4.0 

Clay 

12.5 

36.3 

9.7 

Clay 

16.0 

24.0 

6.8 

Clay 

22.0 

27.0 

6.2 

Clay 

Camouflet 

500 

16.0 

42.0 

9.0 

Clay 

torn 

20.0 

57,0 

14.0 

Clay 

Table  15.  Probable  size  of  shell  craters  in 

compact  soil.  (After  U.S.  Army  1%5.)  Table  16.  Measured  crater  sizes  for  shells  with  delay  fuses 

penetrating  sandy  clay.  (After  U.S.  Army  1965.) 

Swall  Moderate 

penetration  depth  penetration  depth  _ A  verage _  Maximum  recorded 


Shell 

Diameter 

(.ft! 

Depth 

(fo 

Diameter 

(fo 

Depth 

(ft) 

Shell 

Length 

tfo 

Width 

(fo 

Depth 

(fo 

Length 

(.fo 

Width 

(fo 

Depth 

(fo 

75  mm 

4.0 

1,5 

5,0 

3.0 

75  mm 

6 

5 

4 

10 

6 

5 

105  mm 

6.5 

2.5 

7.5 

3.75 

105  mm 

7 

7 

4 

8 

8 

4 

155  mm 

10,0 

4.0 

12.0 

5.0 

155  mm 

12 

11 

6 

14 

12 

8 

8  in. 

11.5 

4.0 

13.5 

5.0 

240  mm 

19 

19 

7 

— 

— 

— 

240  mm 

14.0 

4.0 

15.5 

5.5 

Figure  58.  Dimensions  of  apparent  craters  produced  by  deeply  buried 
bombs.  (Rooke  et  at.  1974.) 
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14  gives  some  values  for  bombs  of  various  sizes 
and  for  various  penetration  depths.  Tables  15  and 
16  give  crater  dimensions  for  shells  of  various  sizes 
with  varying  depths  of  penetration  prior  to  det¬ 
onation.  Trends  of  more  recent  data  were  given  by 
Rookeet  al.  (1974),  as  shown  in  Figure  58.  Craters 
from  deeply  buried  bombs  are  deeper  and  narrow¬ 
er  than  those  from  equivalent  cratering  charges, 
perhaps  because  the  ballistic  penetration  prepares 
a  preferential  venting  path.  Near-surface  detona¬ 
tion  of  shells  apparently  forms  craters  somewhat 
larger  than  those  of  equivalent  stationary  charges, 
probably  because  the  kinetic  energy  at  impact  con¬ 
tributes  to  crater  formation  (inert  projectiles  form 
impact  craters  that  are  quite  similar  to  explosion 
craters). 


EJECTED  FRAGMENTS 
(Flyrock,  ejecta,  missiles) 

In  planning  blasting  operations  or  assessing  ex¬ 
plosion  effects,  it  is  useful  to  know  something 
about  the  range,  size  and  quantity  of  flyrock.  It 
may  also  be  useful  to  have  information  on  impact 
angles,  maximum  heights  reached  by  projectiles, 
and  the  thickness  of  debris  blanketing  the  area 
around  the  crater. 

The  characteristics  of  flyrock  depend  on  the 
charge  weight,  the  charge  depth,  and  the  proper¬ 
ties  of  the  ground  material.  The  range  and  altitude 
of  flyrock  increase  as  charge  size  increases;  the 
total  amount  of  flyrock  and  the  thickness  of  the 
debris  layer  also  increase  with  charge  size.  Loose 
material,  such  as  sand,  obviously  dislodges  more 
easily  than  cohesive  material  like  rock,  but  big 
chunks  of  cohesive  material  can  be  thrown  much 
further  than  small  grains.  A  charge  lying  directly 
on  the  ground  surface  (zero  charge  depth)  acts 
something  like  a  punch,  or  an  impacting  projec¬ 


tile,  kicking  flyrock  up  from  an  annular  zone 
around  the  charge  at  fairly  steep  ejection  angles 
(Fig.  59).  When  the  charge  is  buried,  flyrock  ejec¬ 
tion  angles  tend  to  increase  as  charge  depth  in¬ 
creases,  while  the  horizontal  range  of  projectiles 
decreases.  Eventually,  just  before  charge  depth 
reaches  the  critical  depth  for  camouflet,  most  of 
the  fragments  rise  vertically  and  flop  back  down 
into  the  crater. 

Quantitative  information  on  flyrock  is  sparse 
and  some  of  the  published  relations  are  puzzling. 
In  some  discussions  of  fragment  ejection,  the 
basic  question  of  how  far  flyrock  travels  is  com¬ 
pletely  ignored. 

Simple  tables  of  safe  distance  from  “charges 
placed  on  or  in  the  ground”  have  been  based  on  a 
criterion  of  3(X)  ft/lb (119  m/kg'''’)  for  military 
operations,  with  an  increase  of  the  safe  distance  to 
350  ft/lb for  quarrying  operations  (U.S.  Army 
1967).  For  heavy  surface  charges  (0.5-500  ton), 
Johnson  (1971)  quotes  relations  developed  by 
Vortman  (1967): 

«max  =  ^  "  (23) 

where  is  the  maximum  range  of  ejected  frag¬ 
ments,  W is  charge  weight,  and  /I  is  a  constant  for 
a  particular  charge  shape.  When  is  in  feet  and 
W  is  in  tons  (TNT  equivalent),  A  has  the  value 
1470  for  a  spherical  charge  and  630  for  a  hemi¬ 
spherical  charge  (round  side  up).  The  exponent  0.4 
is  not  explained:  Rooke  et  al.  (1974)  suggest  that 
the  range  of  missiles  from  surface  explosions 
scales  as  ^  (which  is  inconsistent  with  their 
scaling  for  buried  explosions).  For  buried  charges, 
the  maximum  flyrock  travel  is  a  function  of  both 
charge  weight  and  charge  depth.  The  only  pub¬ 
lished  relations  (Fig.  60-61;  see  Vortman  1967, 
Johnson  1971,  Rooke  et  al.  1974)  employ  a  curi¬ 
ous  scaling.  Charge  depth  is  scaled  in  the  usual 


Figure  59.  Trajectories  of  flyrock  for  20-ton  and  100-ton  charges  in  direct  contect  with  the  ground  surface. 


Figure  60.  Scaled  f/yrock  range  as  Figure  61.  Scaled  fly  rock  range  for  charges  in,  and 

a  function  of  scaled  depth  for  above,  hard  rock  and  day  shale.  See  text  for  queries  con- 

charges  buried  in  soil  and  rock.  See  cerning  sixth  rooth  scaling.  (After  Rooke  et  a!.  1974.) 

text  for  queries  concerning  sixth 
root  scaling,  (/if ter  Rooke  et  al. 
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Figure  62.  Maximum  flyrock  range  plotted  against  charge  weight  for  surface  charges 
.set  on  various  ground  materials.  The  shaded  band  indicates  a  cube  root  trend. 
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Figure  63.  Maximum  flyrock  range  plotted  against  charge  depth  in  rock, 
employing  cube  root  scaling  for  both  range  and  depth. 


way  with  respect  to  the  cube  root  of  charge 
weight,  but  flyrock  range  is  scaled  with  This 
implies  a  very  weak  dependence  on  charge  size.  It 
also  conflicts  with  the  finding  that  time  and  veloci¬ 
ty  scale  as  (Cooper  and  Sauer  1977,  Wisotski 
1977),  as  required  by  dimensional  analysis  when 
the  gravitational  acceleration  is  invariant.  Since 
large  fragments  follow  almost  ballistic  trajector¬ 
ies,  elementary  ballistics  gives  range  proportional 
to  the  square  of  velocity,  time,  or  the  product  of 
velocity  and  time,  i.e.  ballistic  range  for  large 
fragments  should  be  proportional  to  or  fV'''’. 

After  examining  the  data  used  by  Vortman 
(1967)  and  later  writers,  it  is  hard  to  accept  the  em¬ 
pirical  scaling  of  flyrock  range  by  a  one-sixth 
power  of  charge  weight.  Only  the  data  for  zero 
charge  depth  (Fig.  62)  can  be  used  to  test  the  hy¬ 
pothesis  directly,  and  these  results  certainly  do  not 
justify  rejection  of  cube  root  scaling.  For  points 
representing  the  same  material  and  charge  type, 
the  exponent  is  about  0.4  (as  in  eq  23).  Cube  root 
scaling  of  all  the  data  for  rock  (Fig.  63)  suggests  a 
plausible  relation  between  maximum  flyrock 
range  and  charge  depth,  with  scatter  similar  to 
that  found  in  typical  data  for  crater  dimensions. 
The  scaled  data  for  soils  (Fig.  64)  show  no  signifi¬ 
cant  correlation  between  flyrock  range  and  charge 
depth,  perhaps  because  these  materials  are  highly 
pulverized  by  the  explosions,  leaving  few  blocky 
projectiles  that  can  be  thrown  to  extreme  range. 

Kinney  and  Graham  (1985),  in  a  discussion  that 


Figure  64.  Scaled  flyrock  range  plotted  against  scaled 
charge  depth  in  soil.  There  is  no  significant  correlation 
for  this  set  of  data  (see  text.) 


seems  to  refer  mainly  to  casing  fragments  from 
surface  explosions,  accept  cube  root  scaling  with¬ 
out  question.  They  give  maximum  missile  range 
^max  as 

=  45  m  (24) 

when  If' is  in  kg  of  TNT  equivalent.  For  safety  dis¬ 
tances  in  bomb  disposal  work,  they  quote 
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Figure  65.  Safety  recommendations  for  bomb  disposal  on  the  ground  surface 
(cube  root  scaling),  together  with  flyrock  range  predictions  for  surface  e.xplo- 
sions  (scaling  exponent  0.4). 


Figure  66.  Variation  of  the  depth  of  deposit¬ 
ed  fragments  with  distance  from  the  explo¬ 
sion.  The  scaling  factor  is  the  cube  root  of 
the  volume  of  the  apparent  crater  (some¬ 
thing  not  likely  to  be  known  before  the 
event).  (After  Andrews  1977.) 


=  120  m/kg'^‘  as  the  desired  distance  and 
=  90  m/kg  '"'  as  the  minimum  distance  (Fig.  65). 

In  a  study  of  flyrock  from  bench  blasts,  Roth 
(1979)  gave  some  consideration  to  flyrock  from 
crater  blasts,  but  did  not  reach  firm  conclusions. 

Some  military  texts  give  more  attention  to  the 
depth  of  ejected  material  on  the  ground  surround¬ 
ing  the  crater  than  to  the  flyrock  range.  Data  are 
given  in  scaled  form,  using  the  cube  root  of  appar¬ 
ent  crater  volume  as  the  normalizing  factor.  Fig¬ 
ure  66  gives  a  relation  between  debris  depth  and 
ground  radius  that  is  based  on  actual  data  points. 
Figure  67  gives  a  relation  extended  to  greater  radii. 
The  uncertainty  in  debris  depth  is  apparently  very 
great  at  the  largest  radii  plotted — more  than  three 
orders  of  magnitude. 

Most  of  the  ejected  material  falls  close  to  the  lip 
of  the  crater.  Figure  68  shows  the  distribution  of 
total  flyrock  volume  with  radial  range  from 
ground  zero.  The  range  is  scaled  with  respect  to 
the  radius  of  the  apparent  crater,  and  it  can  be 
seen  that  for  buried  explosions  about  80%  of  the 
material  is  deposited  within  about  5  crater  radii 
from  ground  zero.  For  near-surface  explosions, 
80%  of  the  flyrock  is  deposited  within  about  7 
crater  radii  from  ground  zero.  Glasstone  et  al. 
(1977)  state  that  the  deposit  of  ejected  material 
from  a  nuclear  explosion  is  continuous  out  to  a  ra¬ 
dius  of  approximately  2  to  3  crater  radii,  and  sug¬ 
gest  that  80%  to  90%  of  the  total  amount  of  eject- 
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Figure  67.  Variation  of  the  depth  of  deposited  frag¬ 
ments  with  distance  from  the  explosion.  Note  the 
very  great  uncertainty  at  the  upper  litnits  of  range. 
fU.S.  .Army  1984.) 


ed  debris  falls  within  this  zone.  The  charge  depth 
is  not  specified. 

The  material  that  travels  farthest  as  flyrock  is 
likely  to  originate  in  a  near-surface  zone  around 
ground  zero.  Conversely,  the  material  that  origin¬ 
ates  close  to  the  boundary  of  the  true  crater  is  like¬ 
ly  to  finish  up  as  fallback.  Figure  69  gives  an  indi¬ 
cation  of  the  likely  travel  distance  for  material 
from  the  crater  of  a  20-ton  explosion  at  a  scaled 
depth  of  approximately  0.8  ft/lb''’  (0.3  m/kg''’)- 
Figure  70  gives  corresponding  information  for  a 
100-ton  explosion  at  a  scaled  depth  of  about  0.2 
ft/lb (0.07  m/kg'^’)- 

For  protection  against  bombardment  by  fly- 
rock,  it  is  useful  to  have  some  information  on 
probable  angles  of  impact.  In  nuclear  weapon 
analyses,  impact  angles  for  ejecta  are  computed, 
making  assumptions  about  ejection  angles  and 
drag  coefficients,  and  using  calculated  impact  vel- 


Figure  68.  Distribution  of  debris  volume  with  distance  from  the 
explosion.  (Rooke  et  al.  1974.) 


Figure  69.  Travel  distance  of  materials  ori¬ 
ginating  in  various  parts  of  the  crater  formed 
by  a  20- ton  explosion.  Material  near  the  top 
of  the  shot  hole  can  be  blown  more  than  25 
crater  radii  from  ground  zero.  Material  that 
originates  near  the  boundary  of  the  true 
crater  is  likely  to  finish  up  as  fallback. 
(Rooke  et  al.  1974.) 
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Figure  70.  Travel  distance  of  material  orif;inat- 
inf;  in  various  parts  of  the  crater  formed  by  a 
lOO-ton  explosion  at  very  shallow  scaled  depth. 
(Rooke  et  al.  1974.) 


Eieciion  Angle  (degrees) 


Figure  71.  Ejection  angles  for  large  fragments 
thrown  out  by  large  surface  explosions.  (U.S.  Army 
1984.) 


ocities  (U.S.  Army  1984).  Ejection  angles  are  as¬ 
sumed  to  range  from  15°  to  75°,  with  the  peak  of 
the  distribution  near  45°  (Fig.  71).  These  compu¬ 
tations  UiC  for  large  surface  explosions.  The  ob¬ 
served  trajectories  for  large  fragments  shown  in 
Figure  59  suggest  that  the  impact  angles  near  max¬ 
imum  range  are  about  50°  to  60°  from  the  hori¬ 
zontal.  There  are  no  available  data  for  buried  ex¬ 
plosions. 

Impact  velocities  have  been  measured  for  frag¬ 
ments  from  surface  explosions,  and  it  has  been 
concluded  that  velocity  scales  with  respect  to  the 
one-sixth  power  of  charge  weight  (Seebaugh  1977, 
Wisotski  1977).  Velocity  varies  with  fragment  size, 
but  the  maximum  fragment  velocities  given 
by  Seebaugh  (1977)  can  be  expressed  as 


F„,a,  =  56.6  fV''*  m/s 
=  \1.2'iW'-  ft/s 


(25) 


with  W  in  tons.  If  this  relation  is  valid,  and  if  the 
particles  it  describes  follow  approximately  ballistic 
trajectories,  then  flyrock  range  should  scale  with 
W'  \  not  fi'  as  shown  in  Figures  60  and  61  or, 
less  seriously,  with  '  as  indicated  by  eq  23. 


GROUND  VIBRATIONS 

Elastic  waves  near  the  ground  surface 

An  earlier  discussion  of  deep  underground  ex¬ 
plosions  dealt  with  the  propagation  and  attenua¬ 
tion  of  elastic  waves.  In  that  discussion,  the  main 
concern  was  with  waves  spreading  spherically  in¬ 
side  solid  ground,  and  with  stress,  strain,  displace¬ 
ment  and  velocity  in  the  ground  material.  How¬ 
ever,  when  blasting  operations  or  weapons  effects 
are  being  considered,  the  main  concern  is  with  dis¬ 
turbance  near  the  surface  of  the  ground,  since  that 
is  where  most  structures  and  human  activities  are 
located.  In  this  section  we  are  concerned  with 
near-surface  regions  that  are  far  enough  from  the 
explosion  for  waves  to  have  attenuated  down  to 
elastic  waves.  In  other  words,  the  “seismic 
region’’  that  is  beyond  the  crater  or  explosion  cav¬ 
ity,  and  also  beyond  the  zone  where  there  is  per¬ 
manent  deformation  of  the  ground. 

The  behavior  of  waves  near  the  ground  surface 
is  very  complicated,  even  for  idealized  conditions. 
An  underground  explosion  first  propagates  body 
waves  in  a  spherical  pattern.  These  body  waves  in¬ 
clude:  1)  the  compression  wave  (also  called 
P-wave,  dilatational  wave,  longitudinal  wave,  or 
irrotational  wave),  and  2)  the  shear  wave  (also 
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called  S-wave,  transverse  wave,  or  rotational 
wave).  As  the  compression  w-ave  passes  through 
an  element  of  ground  material,  it  produces  tran¬ 
sient  compression  and  rarefaction,  with  particle 
motion  in  the  radial  direction  of  propagation.  The 
shear  wave  is  the  transverse,  or  distortional,  com¬ 
ponent  of  the  disturbance,  with  particle  motion 
and  strain  that  is  normal  to  the  radial  propagation 
direction.  The  shear  wave  is  slower  than  the  com¬ 
pression  wave,  roughly  by  a  factor  of  2.  Displace¬ 
ments  in  a  horizontally  propagating  S-wave  are 
resolved  into  horizontal  and  vertical  components 
labeled  SH  and  SV  waves  respectively. 

On  reaching  a  horizontal  ground  surface,  the 
body  waves  reflect  and  they  also  transform  part  of 
the  energy  into  surface  waves,  which  travel  along 
the  surface,  creating  significant  disturbance  to  a 
depth  of  about  one  wavelength.  The  P-wave  and 
the  vertical  component  of  the  S-vvave  couple  and 
create  a  surface  wave  called  the  Rayleigh  wave,  in 
which  particle  motion  follows  an  elliptic  path  in  a 
vertical  plane  that  is  parallel  to  the  radial  propaga¬ 
tion  direction.  The  Rayleigh  wave  is  a  bit  slower 
than  the  shear  wave,  and  its  amplitude  attenuates 
exponentially  with  depth  below  the  surface.  The 
other  type  of  surface  wave  is  a  Love  wave,  in 
which  particle  motion  is  transverse  in  the  horizon¬ 
tal  plane,  i.e.  an  SH  wave.  Love  waves  depend  on 
the  existence  of  one  or  more  layers  parallel  to  the 
surface.  Rayleigh  waves  are  non-dispersive  in  a 
uniform  medium  (i.e.  propagation  velocity  is  inde¬ 
pendent  of  frequency)  but  both  Rayleigh  waves 
and  Love  waves  are  dispersive  in  a  layered  system 
(i.e.  propagation  velocity  varies  with  frequency). 

In  real  situations  the  ground  surface  is  often  ir¬ 
regular  instead  of  being  perfectly  flat,  and  the  sub¬ 
surface  geology  can  give  varying  rock  properties, 
with  multiple  layers  and  interfaces  that  are  not 
necessarily  parallel  to  the  surface.  Wave  propaga¬ 
tion  then  becomes  much  more  complicated  than  it 
is  in  a  uniform  medium  or  in  a  simple  layered  sys¬ 
tem.  Furthermore,  the  explosion  itself  is  not  al¬ 
ways  a  fully  contained  underground  explosion;  it 
can  vent  and  create  air  blast,  so  that  some  ground 
disturbance  is  excited  by  air  blast.  In  commercial 
blasting  where  delays  are  used,  multiple  explo¬ 
sions  occur  in  close  sequence,  sending  out  a  suc¬ 
cession  of  waves  from  different  source  locations. 
Detailed  analysis  of  the  various  waves  and  their 
components  is  important  in  some  areas  of  re¬ 
search,  but  because  of  all  the  complications,  dam¬ 
age  predictions  have  to  be  based  on  a  simplified 
description  of  net  ground  motion. 

In  the  consideration  of  weapons  effects,  the 


problem  of  ground  motion  tends  to  be  more  gen¬ 
eral.  In  military  texts,  the  body-wave  and  shear- 
wave  type  of  disturbance  that  was  just  discussed  is 
called  direct-induced  ground  shock,  since  it  results 
from  direct  coupling  of  energy  to  the  ground.  An¬ 
other  disturbance,  lagging  well  behind  the  direct 
shock,  can  be  produced  by  the  eruption  effects 
during  crater  formation.  This  is  referred  to  as 
crater-induced  ground  shock.  When  explosions 
are  near  the  ground  surface,  or  above  it,  there  is 
another  disturbance  called  airblast-induced ground 
shock.  This  disturbance  is  a  complicated  one, 
since  the  airblast  shock  front  travels  over  the  sur¬ 
face,  giving  a  distributed  source  of  varying  inten¬ 
sity,  with  varying  times  of  origin.  The  following 
notes  cover  only  direct  ground  shock,  with  empha¬ 
sis  on  the  radial  distances  at  which  ground  motion 
is  near  the  limit  for  causing  structural  damage. 

Ground  motion  is  usually  measured  by  some 
type  of  seismograph  set  on  the  ground  surface,  or 
attached  to  the  surface.  Ideally,  it  should  measure 
three  orthogonal  components  of  acceleration 
(preferably  horizontal  radial,  horizontal  trans¬ 
verse,  and  vertical),  and  it  should  record  the  dis¬ 
turbance  as  a  function  of  time.  Some  simple  de¬ 
vices  measure  only  peak  values  of  disturbance. 
The  ground  motion  can  be  described  in  terms  of 
acceleration,  velocity,  and  displacement,  either 
for  a  single  directional  component  or  for  the  re¬ 
sultant  of  the  three  components  (vector  sum). 
When  assessing  the  effects  of  commercial  blasting, 
ground  vibration  is  now  usually  expressed  in  terms 
of  particle  velocity.  For  prediction  of  weapons  ef¬ 
fects,  both  acceleration  and  velocity  are  used.  Fre¬ 
quency  is  a  significant  variable,  but  it  is  not  always 
feasible  to  include  it  in  predictive  schemes.  Meas¬ 
uring  devides  for  monitoring  commercial  blasting 
should  have  a  frequency  range  of  approximately  2 
to  200  Hz. 

Attenuation  of  surface  vibrations 

When  an  explosion  is  buried  and  the  point  of  ob¬ 
servation  is  well  below  the  surface,  wave  propaga¬ 
tion  is  essentially  spherical  and  it  is  reasonable  to 
scale  distance  with  respect  to  the  cube  root  of 
charge  weight  (or  energy  yield),  as  discussed  ear¬ 
lier.  By  contrast,  when  a  concentrated  explosion  is 
underground  but  the  point  of  observation  is  at  the 
ground  surface,  part  of  the  energy  propagates 
spherically  into  the  earth  but  part  spreads  two- 
dimensionally  in  the  surface  layers.  Under  these 
conditions,  neither  cube  root  scaling  (spherical 
spreading)  nor  square  root  scaling  (cylindrical 
spreading)  is  wholly  appropriate.  As  a  first  guess. 
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one  might  expect  cube  root  scaling  to  be  the  best 
approximation  at  small  scaled  radii  from  ground 
zero  for  a  deep  blast,  with  square  root  scaling  giv¬ 
ing  the  best  approximation  at  large  scaled  radii 
from  a  shallow  blast,  where  the  propagation  path 
is  not  far  from  horizontal  and  most  of  the  disturb¬ 
ance  is  caused  by  surface  waves,  which  propagate 
in  two  dimensions. 

Prediction  equations  for  ground  motion  from 
large  buried  explosions  were  given  by  Hughes 
(1968),  Johnson  (1971)  and  Mattes  (1971).  The  at¬ 
tenuation  equations  for  peak  particle  acceleration 
a  and  peak  particle  velocity  v  were  given  in  the 
form 


(26) 

(27) 

where  is  charge  weight,  r  is  radial  range,  and 
the  constants  and  depend  on  the  ground  ma¬ 
terial  at  the  detonation  point  and  at  the  observa¬ 
tion  point.  Rewriting  these  in  terms  of  scaled 
radius: 

a  =  k^{r/fV°^^)~^  (28) 

V  =  (29) 

Thus  the  scaling  is  intermediate  between  cube  root 
( IT*’  ^^)  and  square  root  (fT®  ^),  but  closer  to  cube 
root.  The  relations  are  considered  valid  for  charge 
weights  from  1  to  1000  tons  (TNT  equivalent). 
From  graphical  displays  of  the  equations  (John¬ 
son  1971)  it  appears  that  the  range  of  applicability 
might  be  approximately  7  <  r  <  300  ft/lb using 
cube  root  scaling  for  convenience.  Values  of 
and  ky  for  fully  contained  explosions  are  shown  in 
Table  17. 


Table  17.  Ground  motion  constants  for  large  under¬ 
ground  explosions.  (From  Johnson  1971.) 


Ground  material 
at  detonation 
point 

Ground  material 
al  point  of 
observation 

k  • 

kv* 

High  strength  rock 

Rock 

4,5x10“ 

20  X  10' 

High  strength  rock 

Soil 

9.7x10' 

67  X  10 

Weak  rock 

Rock 

1.8x10' 

8.3  X  10' 

Weak  rock 

Soil 

4.2x10' 

28x10' 

.Soil 

Rock 

0.8x10' 

3.3  X  10' 

.Soil 

Soil 

1.6x10' 

11x10' 

Figure  72.  Attenuation  of  peak  particle 
acceleration  with  scaled  distance  for  var¬ 
ious  combinations  of  ground  conditions. 

See  text  for  notes  on  applicability. 
(Adapted  and  scaled  from  graphs  given 
by  Mattes  1971.) 

Figures  72  and  73  give  graphical  approximations 
of  eq  28  and  29,  using  cube  root  scaling  for  radial 
distance.  These  approximations  should  be  accept¬ 
able  for  charge  weights  in  the  range  0. 1  to  10  tons, 
and  for  fairly  short  distances  from  the  explosion. 

Ground  vibrations  from  commercial  blasting 
operations  are  usually  described  and  regulated  in 
terms  of  peak  particle  velocity,  and  the  main  con¬ 
cern  is  often  vibration  at  considerable  distance 
from  the  blast.  The  U.S.  Bureau  of  Mines  accepts 
square  root  scaling  of  distance  as  the  best  approxi¬ 
mation  under  these  circumstances  (Wiss  and  Line- 
han  1979,  Siskind  et  al.  1980,  Dick  et  al.  1983),  as 
do  other  agencies  (e.g.  Konya  and  Walter  1985). 
The  attenuation  relation  for  peak  particle  velocity 
V  then  takes  the  form 

V  =  kir/W'^f-"  (30) 

where  W  is  the  charge  weight  for  each  separate  ex¬ 
plosion  (charge  weight  per  delay).  In  studies  for 
the  Bureau  of  Mines,  scattered  data  from  a  variety 
of  sources  were  used  to  derive  many  different  val¬ 
ues  of  k  and  n,  and  separate  values  were  obtained 


•  Units  of  and  are  such  that  a  and  v  are  given  in  multiples  of 
g  and  in  in./s  respectively  when  r  is  in  feet  and  W  is  in  tons. 
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(mm/s)  (in./s) 


10°  lo'  10^  10^ 

I  hill _ I  I  I  ililil _ i  .1  -LiliiJ _ I  I  I  jlihl 

10°  lo'  10^  10^  (m/kg'^^) 

Scaled  Distonce 


Figure  73.  Attenuation  of  peak  particle 
velocity  with  scaled  distance  for  various 
combinations  of  ground  conditions.  See 
text  for  notes  on  applicability. 
(Adapted  and  scaled  from  graphs  given 
by  Mattes  1971.) 


Figure  74.  Representative  attenuation  for  particle  velocity  in 
typical  ground  material,  as  determined  from  vibration  studies 
relating  to  commercial  blasting.  (See,  for  example,  Siskind  et 
al.  1980.) 


for  radial,  vertical  and  transverse  components  of 
particle  motion.  No  firm  recommendations  for 
overall  representative  values  of  k  and  n  were  of¬ 
fered,  but  n  was  mostly  in  the  range  1.2  to  1.8. 
The  range  of  data  for  v  was  approximately  0.006 
to  7  in./s  (0.15  to  180  mm/s),  with  {r/W'''^)  in  the 
range  3.5  to  800  ft/lb DuPont  (1980)  accepts  the 
values  n  =  1.6  and  A:  =  160  when  v  is  in  in./s,  r  is 
in  ft,  and  IV  is  in  lb  (per  delay  period  of  8  ms  or 
more).  Konya  and  Walter  (1985)  accept  n  =  1.6 
and  k  -  100  for  the  radial  component  when  the 
same  units  are  used.  The  attenuation  will  obvious¬ 
ly  vary  with  the  type  of  ground  material,  but  Fig¬ 
ure  74  shows  representative  behavior  for  unspeci¬ 
fied  ground  conditions. 
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